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I.  SYNOPSIS  OF  CONTRACT  ACTIVITIES 

It  has  been  a  rather  fruitful  and  eventful  three-year  period  at  the  MIT 
Center  for  Space  Research  with  regard  to  our  ionospheric  and  magnetospneric 
research,  beam-plasma  turbulence  studies,  and  other  related  professional  activi¬ 
ties.  We  have  made  some  definitive  strides  towards  the  basic  understanding  of 
the  relevant  plasma  instabilities,  weak  and  strong  turbulence,  and  wave-particle 
interactions  in  the  Earth's  ionosphere  and  magnetosphere.  The  following  is  a 
summary  of  our  contract  activities: 

1.  Using  the  idea  of  electron-beam  generated  lower  hybrid  turbulence,  we 
were  able  to  demonstrate  how  the  auroral  F-region  ions  can  be  energized  and 
transported  to  the  suprauroral  region  of  the  magnetosphere.  We  were  able  to 
show  that  the  resulting  energetic  ion  conics  and  beams  could  excite  the  observed 
electrostatic  i on  cyclotron  waves  and  ion  cyclotron  harmonic  waves  at  hiqher 
altitudes.  We  have  published  four  papers  and  presented  several  invited  lectures 
and  special  topics  seminars  on  this  subject. 

Our  research  has  been  aided  by  the  useful  interactions  with  the  able  plasma 
theorists.  Dr.  J.  R.  Jasperse  at  the  Air  Force  Geophysics  Laboratory,  Drs.  B.  Basu 
and  J.  Retterer  of  the  Space  Data  Analysis  Laboratory,  Boston  College,  and 
Drs.  R.  L.  Lysak  and  M.  K.  Hudson  of  the  University  of  California  at  Berkeley. 
Dialogues  with  the  experimental  groups  (Drs.  J.  D.  Winningham  and  J.  Burch  at 
the  Southwest  Research  Institute,  Dr.  D.  Klumpar  of  the  University  of  Texas  at 
Dallas,  Dr.  P.  Kintner  of  the  Cornell  University,  Dr.  J.  Sharber  of  the  Florida 
Institute  of  Technology,  and  Dr.  M.  Temerin  at  the  University  of  California  at 
Berkeley)  responsible  for  the  particle  and  wave  measurements  of  the  DE,  ISIS-2, 
and  S3-3  satellites  and  other  high  altitude  sounding  rockets  were  instrumental 
in  the  development  of  our  basic  theoretical  model. 
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2.  We  have  made  some  relevant  plasma  wave  calculations  related  to  the 
interesting  beam-plasma-discharge  (BPD)  phenomenon  produced  by  an  electron  beam 
iniected  in  a  vacuum  rhairhp'''  under  amMart  niac.nv  m nnoUr  •<  ^nr^U'-rc 
simulating  the  ionospheric  environment.  Similar  ideas  on  the  plasma  waves  and 
the  associated  wave-particle  interactions  induced  by  artificially  produced  charged 
(ion  and  electron)  particle  beams  injected  in  the  ionosphere  aboard  sounding 
rockets  have  also  been  developed.  We  have  enjoyed  the  discussions  with  Mr.  H.  Cohen 
and  Dr.  W.  Denig  at  the  Air  Force  Geophysics  Laboratory,  Dr.  P.  Kellogg  at  the 
University  of  Minnesota,  Dr.  P.  Kintner  at  the  Cornell  University,  and  Drs.  R.  Arnoldy 
and  T.  Moore  of  the  University  of  New  Hampshire  at  Durham,  in  connection  with  the 
various  experimental  facets  of  these  fascinating  beam-generated  plasma  phenomena 

in  space  and  in  the  laboratory. 

3.  We  have  made  exact  nonlinear  plasma  wave  calculations  capable  of  describ¬ 
ing  sawtootn  and  double-spike  shaped  finite  amplitude  electrostatic  ion  cyclotron 
waves  similar  to  those  observed  in  the  suprauroral  region.  We  have  also  demonstrated 
the  existence  of  a  class  of  nonlinear  electrostatic  ion  cyclotron  cavitons.  The 
structures  of  these  cavitons  bear  striking  resemblance  of  the  paired  electrostatic 
shocks  frequently  detected  near  1  above  the  auroral  zones.  We  have  Dublisbed 

two  papers  and  presented  several  seminars  on  this  subject. 

4.  Jointly  with  Dr.  J.  R.  Jasperse  at  the  Air  Force  Geophysics  Laboratory, 
we  have  been  studying  the  general  properties  of  the  collision  operators  that 

ar«  relevant  to  the  Bol tzmann-Fokker-Pl anck  equation  for  ionospheric  applications 
We  nave  presented  several  papers  on  this  subject  at  the  AGU  and  APS- Di vision  of 
■/sics  Meet  Jigs. 

5  Dr.  H.  Koons  of  the  Aerospace  Corporation  and  Mr.  H.  Cohen  at  the  Air 
'>onhvsics  Laboratory  have  provided  us  with  some  very  in*<  re'  *  in  na«ma 
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wave  data  produced  by  artificially  injected  electron  and  ion  beams  carried  aboard 
the  Scatha  satellite  which  has  geostationary  orbits.  One  of  the  surprising 
phenomena  indicated  by  the  data  is  the  presence  of  electromagnetic  emissions 
at  roughly  3/2  times  the  electron  gyrofrequency.  We  have  found  a  direct  non¬ 
resonant  electron-beam  generation  mechanism  which  might  be  responsible  for  such 
plasma  wave  emissions.  This  work  was  carried  out  in  collaboration  with 
Ms.  D.  Donatelli  of  the  Physics  Department,  Boston  College,  and  Dr.  W.  Burke 
of  the  Air  Force  Geophysics  Laboratory. 

6.  Throughout  the  Earth’s  magnetosphere  and  ionosphere  (as  well  as  near 
the  bow  shock  and  in  the  solar  wind),  low  frequency  electrostatic  noises  have 
been  detected.  A  number  of  these  emissions  had  been  identified  as  ion  acoustic 
waves.  We  found  that  when  the  ion  temperature  is  high,  the  plasma  medium  generally 
cannot  support  ion  acoustic  waves.  Some  of  these  emissions  can  be  identified  as 
lower  hybrid  waves.  Sometimes,  the  electric  field  noises  are  accompanied  by 

low  frequency  electromagnetic  fluctuations.  Such  emissions  generally  propagate 
almost  normal  to  the  ambient  magnetic  field  lines.  We  found  that  they  are,  in 
fact,  electromagnetic  lower  hybrid  waves  or  hybrid  whistlers.  This  work  is  done 
in  collaboration  with  Dr.  E.  Marsch  of  the  Max-Planck  Institut  fiir  Aeronomie. 

We  have  published  two  papers  in  the  Geophysical  Research  Letters  and  the  Journal 
of  Geophysics  Research  and  presented  a  number  of  invited  talks  on  this  subject. 

7.  Using  the  path-integral  and  renormal ization-group  approach,  we  have 
formulated  a  procedure  for  calculating  the  correlation  and  response  functions 

of  stochastic,  nonlinear,  nonequilibrium  systems.  These  ideas  should  have  direct 
applications  to  the  analytical  description  of  strong  plasma  turbulence.  We  are 
currently  applying  these  ideas  to  the  lower  hybrid  turbulence  and  stochastic 
plasma  heating.  Dr.  D.  Vvedensky  of  the  Imperial  College  of  London  and 
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Dr.  J.  F.  Ni col  1  of  the  University  of  Maryland  joined  us  in  developing  the  above 
ideas.  We  have  published  two  papers  in  the  Physics  Letters  A,  and  Physical 
Review  A  on  this  subi<=rT 

8.  Magnetic  reconnections  are  believed  to  take  place  at  the  magnetotail 
and  on  the  dayside  of  the  magnetosphere.  It  has  been  suggested  that  the  magnetic 
reconnection  process  might  be  related  to  the  onset  of  magnetic  substorms.  The 
particle  distributions  in  the  magnetosphere  is  generally  believed  to  be  anisotropic. 
We  have  studied  the  anisotropy-dri ven  magnetic  reconnection  phenomenon.  Two 
papers  on  magnetic  reconnection  and  related  instabilities  have  been  published 

in  the  physics  journals.  This  work  is  done  in  collaboration  with  Dr.  B.  Basu 
of  the  Space  Data  Analysis  Laboratory  at  the  Boston  College  and  Dr.  T.  Antensen 
of  the  University  of  Maryland. 

9.  During  the  contract  period,  we  organized  or  participated  in  organizing 
three  national  and  international  conferences: 

1.  AGU  (American  Geophysical  Union)  Chapman  Conference  on  "Waves 
and  Instabilities  in  Space  Plasmas,"  held  in  Denver,  Colorado 
on  August  7-9,  1979. 

2.  "Plasma  Astrophysics"  sponsored  by  the  International  School 
of  Plasma  Physics  and  held  at  Varenna  (Como),  Italy  on 
August  27  -  September  8,  1981. 

3.  The  1982  "International  Conference  on  Plasma  Physics"  held  in 
Goteborg,  Sweden  on  June  9-15,  1982. 

In  January  each  year,  in  cooperation  with  Dr.  J.  Jasperse  of  the  Air  Force 
Geophysics  Laboratory,  we  have  organized  a  one-day  symposium  on  the  "Physics 
o+  .pace  Plasmas"  as  one  of  MIT's  annual  IAP  (Independent  Activity  Period) 

These  symposia  were  sponsored  by  the  Center  for  Spare  Pr search  <t  ’ 
arid'rts -Director,  Dr.  H.  S.  Bridge.  The  topics  discussed  at  the  symposia  covered: 

11 i  It.petary  ionospheres- and  magnetospheres,  the  solar  wind,  wave-ovt  irle  in  cr 
1  i 'id  particle  acceleration  mechanisms ,  aurora  rlasma  r>'ivsi.--  .  am  ..-.ip.  »  •  < 


reconnection  processes  in  space.  The  proceedings  of  these  symposia  are  published 
as  hard-covered  books  by  the  Scientific  Publishers,  Inc.,  Cambridge,  Massachusetts. 
During  the  three-year  contract  period,  we  chaired  several  space  plasma  sessions 


and  gave  a  number  of  invited  talks  at  the  annual  AGU  and  APS  (American  Physical 
Society)  meetings  as  well  as  other  national  and  international  conferences.  We 
have  presented  numerous  seminars  on  our  research  findings  at  various  major 
universities  and  scientific  institutions. 

The  content  of  this  report  is  oraanized  as  follows.  In  Sect.  a  II-VIII, 
we  discuss  in  detail  the  major  contributions  of  our  research  acti\  <s.  In 
Section  IX,  we  list  the  titles  of  the  published  technical  and  scii  *-'c  papers. 
In  Section  X,  we  list  the  titles  of  the  papers  presented  at  scient,,.c  meetings 
and  conferences . 
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II.  ION  CONICS  IN  THE  TOPSIDE  IONOSPHERE  AND  THE  SUPRAURORAL  REGION  OF  THE 

MAGNETOSPHERE  _ _ _ _ _ _ _ 

A.  INTRODUCTION 

Recent  satellite  and  rocket  data  on  particle  distributions  collected  from 
the  topside  ionosphere  and  the  suprauroral  region  indicate  that  ions  (H+  and  0+) 
of  ionospheric  origin  are  accelerated  transverse  to  the  geomagnetic  field  lines 
(Sharp  et  al.,  1977;  Lysak  et  al . ,  1980;  Klumpar,  1979;  Ghielmetti  et  al . ,  1978; 
Gorney  et  al.,  1981).  The  energy  range  of  the  upward  flowing  ions  rarely  extend 
beyond  500  eV  at  lower  altitudes  (<5,000  km)  and  the  distributions  are  usually 
"shallow  conics"  with  pitch  angles  average  between  90°  to  140°.  (Without  loss  of 
generality,  we  restrict  our  discussions  to  the  northern  hemisphere.)  At  altitudes 
above  5,000  km,  the  ions  can  be  accelerated  beyond  the  kilovolt  range  (Kintner 
et  al.,  1979)  and  the  average  pitch  angles  become  closer  to  180°. 

Throughout  the  suprauroral  region,  ac  electric  field  measurements  have 
detected  plasma  waves  (Gurnett  and  Frank,  1972;  Lae'pere  and  Hoffman,  1976; 

Mozer  et  al.,  1979;  James,  1976)  dominated  generally  by  lower  hybrid  (LH)  modes. 

At  altitudes  above  5,000  km,  where  kilovolt  ion  conics  and  beams  have  been 
observed,  electrostatic  ion  cyclotron  modes  are  also  detected  (Kintner  et  al., 
1979;  Temerin  et  al.,  1979). 

It  is  known  that  energetic  electrons,  injected  into  the  suprauroral  reoion, 
from  the  plasma  sheet,  are  accelerated  by  field-aligned  dc  electric  fields  or 
double  layers  towards  the  ionosphere.  The  energies  of  these  electron  beams  can 
reach  the  keV  range.  A  sequence  of  electrostatic  modes  can  be  excited  by  wave- 
particle  resonances  involving  the  positive  slope  region  of  the  electron  distribu¬ 
tions  (Coppi  et  al . ,  1974;  Maags,1976;  Maggs  and  Lotko,  1981;  Papadopoulos  and 
Palmadesso,  1976).  The  modes  in  the  frequency  range  of  the  ion  plasma  frequency 
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(LH  modes)  can  resonate  simul taneously  with  both  the  electron  and  ion  popula¬ 
tions  (Coppi  et  al.,  1976),  thus  transferring  energy  from  the  electrons  to  the 
i  ons . 

This  ion  acceleration  process  is  particularly  efficient  at  altitudes  between 
1,000  Km  and  5,000  km  where  LH  waves  are  generally  broad-band  and  the  wave  intensity 
is  relatively  high  (Chang  and  Coppi,  1981).  It  can  be  shown  that  ions  (H+  and 
0+)  of  ionospheric  origin  can  be  accelerated  upwards  along  the  field  lines  and 
form  ion  conics . 

6 .  EXCITATION  OF  ELECTROSTATIC  LOWER  HYBRID  WAVES  BV  ENERGETIC  ELECTRON  BEAMS 

We  consider  electrostatic  modes  with  E  =  -  7<J>  and  T  =  exp(-iu>t  +  ik„Xn  4- 

ik^xj,  where  the  subscripts  >•  and  x  indicate  the  parallel  and  perpendicular 

components  of  the  wave  vector  k  to  the  geomagnetic  field.  In  addition, 

-2  2  -2 

■■  "<•  '  and  <<  k  «  ?_  where  9..  -  eB/in-c  and  0  =  eB/ni  c  represent 

1  c  1  X  6  I  1  c  t5 

the  ion  and  electron  qyrofrequencies ,  while  p..  =  v.j/f;.  and  pg  =  ve/fi  represent 
the  average  ion  and  electron  gyroradii.  Thus,  we  need  only  consider  the  electron 
dynamics  along  the  magnetic  field  line  while  the  ions  can  be  treated  as  unmaqnet i zed. 
Consequently,  we  consider  the  following  one-dimensional  electron  and  icn 

mole1  distributions: 

F  =  n  (m  / 2ttT  )1/2  exp  (-v2/v  2) 
e"  eo  e  eo  v  '  eo 

+  n  u(fiU/2mT  ,)^2  exp  [  -  ( v  i.  -u  ,  ) 2  /  v  fl.  Ma' 
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where  u  .  represents  the  electron  beam  velocity,  v  -  2T  /m  ,  v  f  =  2T  ,/m  , 
eo  '  eo  e„o  e  eb  e»b  e 

2  2 

vio  =  ^T^/m^ ,  vih  =  2Tiih/mi,  vA  =  kA  •  y/kA,  the  subscript  "o"  refers  to  the 
main  body  of  both  the  electron  and  ion  distributions,  while  the  subscript  "b" 
refers  to  the  beam  electrons,  and  *he  subscript  "h"  to  the  hot  ion  component  of 
the  ions  produced  by  the  excitation  of  the  considered  LH  modes. 


For  u)/k„veo  >>  1,  and  w/kv^  »  1,  the  dispersion  relation  becomes: 
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Here,  we  have  assumed  <<  wcg,  oi  .  >>  u>ci.  which  are  conditions  usually 
satisfied  in  the  suprauroral  region.  In  (2),  we  have  defined: 


For  (neb’  nih}  <<  neo  ~  nio’  we  have; 
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with  oj  =  +  iy,  where  ujp..  =  (47in^oe  /m^ )  is  the  ion  plasma  frequency.  For 

co/k,,  ?  u^  and  u/kA  >  v^,  these  waves  can  resonate  simultaneously  with  the 

electrons  and  ions,  thus  transferring  energy  and  momentum  from  the  beam  electrons 

to  the  hot  ion  component  (Coppi  et  al.,  1976).  Thus,  from  (4)  we  conclude  that 

u>r  of  the  LH  modes  in  the  suprauroral  region  is  of  the  order  of  and  it  is 

cut  off  from  below  at  w  .. 

Pi 

We  can  argue  that  the  considered  instability  becomes  quasi-stationary  when 
the  characteristic  parameters  of  the  beam  electron  distribution  and  of  the  heated 
ions  make  the  expression  for  y  obtained  from  the  linearized  stability  theory 
vanish.  Then, 
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As  indicated  by  the  laboratory  experiments  tnat  have  led  to  identify  a  "slide- 
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that  this  process  can  provide  a  very  efficient  mechanism  for  producing  a  hot 

2  2 

component  in  the  ion  distribution.  Since  k„  <<  kx,  and  the  relevant  resonant 
condition  is  w  -  kxvx) .  =  o,  the  ions  receive  energy  from  the  excited  mode  in 
the  transverse  direction  to  the  magnetic  field  lines. 


C.  QUAS I L I  NEAR  ION  EVOLUTION  EQUATION 

Because  the  lower  hybrid  modes  observed  in  the  lower  suprauroral  region  are 
generally  broad-band,  we  can  use  the  quasilinear  expression  to  describe  the  ion 
heating  process.  For  unmagnetized  ions,  the  quasilinear  wave-particle  interaction 
term  is: 


(Sfi  g  3^ 

6t~ ' QL  =  3y  ‘ 
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is  the  diffusion  coefficient.  In  this  expression,  f.  is  the  ion  distribution 


function,  and  S  is  the  spectral  density  of  the  electrostatic  potential  fluctuations. 


To  evaluate  the  quasilinear  term,  we  need  to  know  the  spectral  density  of 
tie  •(! r-hyhrid  turbulence.  A  self-consistent  calcic  at  ion  of  S  ,  ( s  • 1 '  is  com- 
i'l  >  by  the  problem  of  understanding  the  saturation  mechanism  of  the  lower 

’  .•!  :I  i-'  /.  whether  it  occurs  by  craved:  •  i  '  -.j  r<  •: . 

’‘'■‘Step  that  problem,  we  w i 1 1  use  experimental  1  v  observed  amplitudes  of  the 
■id  waves .  Although  low®r  h/hvid  f  u*"hi»l  ;iv»  >'s  to  '  •  i->v 


in  frequency,  little  is  known  concerning  the  wave-vector  spectrum.  A  variety 


of  different  assumptions  for  S  (k,w),  however,  all  lead  to  estimates  of  the 

y 

diffusion  coefficient  which  are  similar  in  order  of  magnitude.  To  this  level 
of  accuracy,  we  have: 


'm.  E 
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where  is  the  power  spectrum  of  the  electric  field  (in  units  of,  e.g., 

2 

(v/m)  /Hz).  Because  the  frequency  spectrum  typically  extends  from  to  .  to  a 

2 

few  wpi.,  we  take  Sg.  z  |E|  /upi,  where  J E |  is  the  observed  wave  field  amplitude. 

To  calculate  DA  as  a  function  of  velocity,  we  would  again  need  the  unobservable 
wave-vector  spectrum  of  the  turbulence;  for  the  sake  of  simplicity,  we  will  assume 
that  D,  is  constant  above  a  minimum  resonant  velocity  va  .  and  that  D,  is  zero 
below  this  velocity. 

In  addition  to  the  wave-particle  interaction  described  above,  the  geomagnetic 
field  and  the  DC  electric  field  also  play  important  roles  in  the  evolution  of  the 
ions  in  the  suprauroral  region.  Thus,  with  the  approximate  quasilinear  term,  the 
kinetic  equation  for  the  ion  distribution  function  is: 
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where  s  denotes  position  along  a  field  line.  We  wish  now  to  study  the  evolution 
of  the  ions  by  solving  this  equation. 


D.  SAMPLE  ION  CALCULATIONS 

A  Monte  Carlo  analysis  of  the  spatially  dependent  equation  (10)  of  the 
ion  distribution  will  be  considered  in  Sec.  HE.  Here,  we  consider  the  evolution 
of  the  parallel  and  perpendicular  energies  of  sample  ions.  In  the  absence  of  the 
LH  modes,  the  ions  could  have  been  tracked  along  a  given  field  line  according  to 
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their  unperturbed  orbits.  On  the  other  hand,  the  effect  of  the  LH  modes  produces 
anomalous  rates  of  increase  of  the  ion  perpendicular  energy  ^  and  parallel 

momentum  P._).  i  t.  These  can  be  estimated  usinq  rhp  nuasilinear  diffusion  *onr., 
6f./<St)  .  Tnus,  we  have 
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where  w^is  the  ion  perpendicular  energy,  v^  is  the  ion  parallel  velocity,  and 
E„  is  the  field-aligned  dc  electric  field. 

Let  us  consider  first  the  evolution  process  of  an  ion  in  the  absence  of 
E„.  Because  of  the  resonance  condition  for  excitation  of  the  LH  modes  by  the 
downward  flowing  electron  beam,  the  direction  of  the  parallel  component  of  the 
phase  velocity  of  the  LH  modes  at  generation  is  downward.  Initially,  when  the 
value  of  is  relatively  small,  the  wave  drag  P-j„)turb  can  be  larger  tnan  the 
mirror  magnetic  field  term  in  (lib).  Since  the  net  momentum  transfer  rate  is 
negative  for  this  case,  the  ion  will  first  acquire  a  small  downward  velocity 
(if  we  assume  the  initial  parallel  velocity  is  zero).  At  this  time,  the  perpen¬ 
dicular  energy  is  much  larger  than  the  parallel  component.  This  means  that  the 
nitch  angle  of  the  ion  is  close  to  90°.  Eventually,  the  perpendicular  energy 
is  increased  to  a  sufficiently  high  magnitude  by  the  LH  modes  for  the  mirror 
for'  •••  n  the  parallel  equation  of  motion  to  overcome  the  wave  drao  term.  The 
i'  r  in i  t  velocity.  will  then  begin  to  decrease  in  magnitude  until  the  motion  cf 
‘  • re  reused. 

h  rind  from  (11a,  b)  that  for  a  root  mean  square  electric  field 
■  1  mV/m  and  typical  densities  in  I  hi*  lower  ‘•uwaurora  I  ren  .  ’>0,  >, 
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H+-  and  0+-ions  can  be  accelerated  from  1  eV  to  tens  or  1 00  eV  in  an  altitudinal 
change  of  the  order  of  100  km.  After  the  ions  are  energized  above  100  eV,  they 
begin  to  gather  upward  momentum  and  the  mirror  conversion  term  in  (11a)  can  no 
longer  be  neglected.  However,  the  drag  term  is  now  much  smaller  than  the  magnetic 
mirror  force  and  can  be  neglected  in  (11b).  Generally,  the  ions  can  be  energized 
to  several  hundred  eV  before  they  are  convected  beyond  the  5,000  km  altitude  and 
the  entire  evolution  process  lasts  for  a  period  between  30-40  seconds  to  four  or 
five  minutes.  (An  illustrative  example  is  given  in  Fig.  1.)  Of  course,  this 
acceleration  process  depends  on  the  existence  of  LH  waves.  These  waves  may 
exist  in  an  altitudinal  range  of  tens  to  hundreds  of  km  or  even  thousands  of  km. 

By  following  the  trajectories  of  the  sample  ions  with  different  initial 
conditions,  we  can  construct  graphically  the  accelerated  ion  distribution  as 
it  evolves  upwards.  Generally,  itattainsa  butterfly-like  configuration  (Fig.  2) 
resembling  those  measured  experimentally  (Mizera  et  al . ,  1981). 

When  the  ions  propagate  across  a  kilovolt  electrostatic  shock  or  a  series 
of  double  layers,  they  tend  to  become  field-aligned.  The  combination  of  such 
keV  ion  distributions  and  of  that  of  the  background  ions  can  lead  to  the  excita¬ 
tion  of  electrostatic  ion  cyclotron  modes  (Kintner,  1980).  Thus,  the  LH  accelera¬ 
tion  mechanism  may  provide  a  possible  explanation  of  the  simultaneous  observation 
of  the  electrostatic  ion  cyclotron  modes  with  the  keV  ion  distributions  in  the 
suprauroral  region,  particularly  at  altitudes  above  5,000  km  (Kintner  et  al . ,  1979). 

E.  A  MONTE  CARLO  TECHNIQUE  FOR  THE  ION  KINETIC  EQUATION 

Because  of  the  inherent  importance  of  its  spatial  structure,  the  time- 
dependent  solution  of  the  kinetic  equation  for  the  ion  distribution  function 
(10)  is  fully  a  four-dimensional  problem--too  compl icated  to  be  solved  by  standard 


finite-difference  techniques.  Instead,  we  adopt  a  model  akin  to  a  particle 
simulation,  in  which  the  effects  of  wave-particle  interactions  are  described 
bv  us  inn  a  Montp  Carl'"'  tpchnioue .  Because  +  hp  ir>n<-  whirh  rpconan+lv  infract 
with  the  lower  hybrid  waves  are  on  the  tail  of  the  distribution  function,  they 
are  few  in  number.  Unless  they  excite  instabilities,  their  effect  on  the  rest 
of  the  particles  is  therefore  small.  Thus,  we  need  not  create  a  completely 
self-consistent  model  of  the  plasma;  instead,  we  can  impose  fields  and  turbulence 
and  study  how  these  affect  the  resonant  ions. 

We  begin  with  a  segment  of  field  line,  from  height  s  =  sL  to  s^.  Alonq 
this  segment,  we  specify  the  geomagnetic  field  strength  B(s),  the  dc  electro¬ 
static  potential  <t>(s),  and  the  magnitude  of  the  quasil  inear  diffusion  coefficient 
DiCs).  The  distribution  of  the  resonant  ions  is  approximated  by  using  a  large 
number  of  individual  particles.  Their  velocity  distribution  is  initially  a 
Maxwellian.  (Actually,  because  ions  with  velocities  below  v^  min  will  not 
interact  with  the  lower  hybrid  turbulence,  we  do  not  bother  to  load  them  into 
our  simulation).  The  particle  density  is  assumed  to  be  a  power-law  function 
of  altitude;  the  exponent  -2  gives  good  agreement  with  Maeda's  (1975)  model  of 
plasma  density  in  the  suprauroral  region. 

The  calculation  proceeds  by  following  the  orbits  of  the  particles  with 
time.  Because  we  are  interested  in  changes  occurring  on  scales  much  larger 
than  tne  size  of  a  gyro  orbit,  we  need  not  integrate  the  equations  of  motion 
ir  f‘.<  tail.  Instead,  we  follow  the  motion  of  the  particle  guiding  .enters  alonn 


If  for  a  moment,  we  neglect  the  effect  of  wave-particle  interactions,  the  velocity 
of  a  particle  as  a  function  of  position  can  be  found  by  assuming  the  invariance 
of  the  particle's  energy  and  of  its  first  adiabatic  invariant  over  a  short  time 
step.  In  this  way,  the  equation  of  motion  could  be  integrated  step  by  step. 

When  we  include  wave-particle  interactions,  we  need  to  modify  this  scheme 
slightly.  In  addition  to  the  change  in  velocity  caused  by  the  particle's  motion 
in  the  systematic  fields,  a  change  in  velocity  due  to  exchange  of  energy  with 
the  turbulent  fields  must  also  be  included  during  each  time  step.  Because  the 
interaction  is  modeled  as  a  diffusion  process,  the  first  approximation  to  the 
probability  of  a  given  change  in  velocity  during  an  interval  At  is--by  the 
central  limit  theorem--a  Gaussian  function  of  Av  in  one  dimension.  The  variance 
of  this  Gaussian  is  2DxAt;  because  there  is  no  drag  term  in  the  kinetic  equation, 
its  mean  is  zero.  In  the  present  case,  diffusion  occurs  in  two  dimensions,  in 
the  plane  perpendicular  to  the  magnetic  field.  Therefore,  the  increment  in  the 
square  of  the  transverse  velocity  is  given  by: 


A(vx^)  =  2vx  A  cos®  +  A^, 


where®is  a  random  angle,  o  <®<  2tt,  and  A  =  (Av  +  Av  )  ,  where  Av  and 

x  y  x 

AVy  are  both  Gaussian  random  variables  with  variance  2DxAt.  Such  a  random  incre¬ 
ment  is  added  to  the  velocity  of  each  particle  following  each  time  step. 


F.  THE  FORMATION  OF  ION  CONICS 

It  is  expected  that  transverse  heating  of  ions  by  lower  hybrid  modes  at 
low  altitudes  will  lead  to  the  creation  of  conic  pitch  angle  distributions  seen 
at  higher  altitudes  in  the  suprauroral  region  (Chang  and  Coppi,  1981).  The 
increase  in  perpendicular  energy  first  creates  a  narrow  pitch-angle  distribution 
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peaked  at  90°.  Then,  as  the  particles  mirror  and  move  up  the  field  line,  the 
variation  of  the  geomagnetic  field  will  cause  some  of  the  perpendicular  energy 
to  be  transferred  to  parallel  motion,  systematically  increasing  the  pitch  anqles 
and  creating  the  conic  distribution. 

In  the  calculations  presented  here,  the  following  parameters  were  adopted. 

Over  the  altitude  range  from  1,000  km  to  5,000  km,  we  followed  H+  and  0+  ions, 
with  a  primordial  temperature  of  1  eV.  No  field-aligned  potential  drop  was  included. 
We  assumed  that  the  wave  amplitude  at  1,000  km  was  50  mV/m,  consistent  with  S3-3 
observations  (Temerin  et  al.,  1981).  The  ion  plasma  frequency,  necessary  to 

3 

calculate  Dx,  was  taken  to  be  4tt  x  10  rad/s;  the  minimum  resonant  velocity 
was  the  thermal  velocity.  We  considered  two  different  cases  for  the  altitudinal 
range  of  the  lower  hybrid  turbulence:  (1)  narrowly  distributed  turbulence,  with 
Dx  nonzero  only  in  the  lowest  500  km  of  the  simulation,  and  (2)  widely  distributed 
turbulence  with  Dx  constant  over  i ,000  km.  The  actual  spatial  distribution  is 
poorly  known,  although  evidence  shows  (Temerin  et  al.,  1981)  that  broad-band 
lower  hybrid  waves  will  be  found  below  the  electrostatic  shock  or  in  the  double 
'ayer  region  whenever  electron  beams  are  detected. 

As  the  simulation  progresses,  we  allow  a  steady  state  to  be  established 
:>y  replacing  every  particle  which  leaves  the  simulation  by  a  particle  picked 
at  random  from  the  primordial  distribution.  We  will  illustrate  this  with  phase 
space  snapshots  taken  after  this  steady  state  has  been  established. 

Let  us  consider  the  model  of  heating  with  a  na-row  range  r-  ‘urbulence 
M",’  In  this  case,  the  heating  and  folding  processes  work  almost  independent iy 
•  .  *her  The  steady  state  has  been  » •earned  t  **<••  -»h  •  t  .  •  i-  ‘  *•;. 

ing  t y  rhe  turbulence.  We  find  that  the  resonant,  ions  have  been  heated  to  a 
ii  *-ist.ic  energy  of  a  few  hundred  eV.  Their  pnch-  ar*ole  •Jirt  •  r  i  c. 
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sharply  collimated  at  all  altitudes,  as  it  folds  from  90°  at  1,000  km  to  about 
150°  at  5,000  km. 

In  contrast,  when  diffusion  occurs  over  a  wide  range  of  altitude,  the  H+-ion 
conic  structures  are  more  diffuse.  The  continued  transverse  heating  destroys  the 
collimation.  Because  of  the  wider  extent  of  the  heating  region,  particles  are 
also  accelerated  to  higher  energies,  typically  of  order  1  keV  for  H+-ions  in 
this  case  (Fig.  3).  The  increase  in  energy  is  not  greater  because  the  heating 
is  self-limiting:  as  a  particle  gains  energy,  it  moves  out  of  the  heating  region 
more  quickly.  The  broader  conics  calculated  in  the  second  model  are  much  like 
the  conics  seen  by  S3-3  (Mizera  et  al.,  1981).  Similar  results  are  obtained 
for  0+-ions.  Even  though  the  heating  rate  for  a  given  wave  spectrum  is  lower 
because  of  the  larger  mass,  0+-ions  can  be  heated  to  comparable  energies  (as 
for  H+-ions).  The  reason  is  that,  because  they  are  heavier,  0+-ions  stay  in 
the  heating  range  much  longer  than  H+-ions. 
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Illustrative  example  of  the  energy  and  pitch  angle  evolution  of  an  H  -ion 
in  the  lower  suprauroral  region.  SE  =  10-6  (V/m)2/Hz  and  W0  =  1  eV  at  the 
initial  altitude  h0  =  500  km.  Sf  is  assumed  to  decrease  at  the  rate  of  5°f 
per  500  km  increase  of  altitude. 


Schematic  contours  of  the  accelerated  port 
function. 
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Monte  Carlo  calculation  of  the  evolution  of  H  -ions.  Acceleration  ranoe 
1  ,000  -  2,000  km  (S[r  =  12.5  x  10"?  (V/m)^/Hz).  Elapsed  time:  50  sec. 

A1 titude:  3,000  km. 
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III.  EXACT  NONLINEAR  SUPERSONIC  PERIODIC  ELECTROSTATIC  ION  CYCLOTRON  WAVES 
OF  ARBITRARY  AMPLITUDE _ 

A.  INTRODUCTION 

Spiky  and  sawtooth  shaped  nonlinear  periodic  electrostatic  ion  cyclotron 
(EIC)  waveforms  have  been  detected  near  1  R^  in  the  suprauroral  region  of  the 
magnetosphere  via  electric  field  measurements  by  the  S3-3  satellite  (Temerin 
et  al.,  1979).  Approximate  analysis  (Chaturvedi,  1976)  and  direct  numerical 
computations  (Temerin  et  al.,  1979)  have  shown  that  such  waveforms  are  admissible 
traveling  wave  solutions  to  the  nonlinear  EIC  wave  equations. 

We  present  here  (Chang,  1982)  a  general  analytical  solution  of  nonlinear 
EIC  waves  of  arbitrary  amplitude.  Exact  periodic  waveforms  are  obtained  for 
supersonic  and  hypersonic  propagations. 

B.  BASIC  EQUATIONS 

We  assume:  (i)  quasineutrality,  n.  *  n  =  n,  (ii)  cold  ions,  T.  «  T  , 

16  16 

(iii)  Boltzmann  electron  distribution,  ng  =  nQ  exp  ( ed>/  kTg ) ,  and  (iv)  a  constant 
applied  uniform  magnetic  field,  BQ.  The  electrostatic  ion  cyclotron  waves  proDa- 
aate  essentially  perpendicular  to  the  field  lines.  Since  the  perpendicular 
electron  motion  and  the  parallel  ion  motion  are  not  important  for  electrostatic 
ion  cyclotron  propagations,  we  need  to  consider  only  the  two  perpendicular 
equations  of  motion  for  the  ions: 

3u/3t  +  17  3u/3x  =  -(e/m.j ) 30/3x  +  eB  v/m..c,  (la) 

3v/3t  +  IT  3v/3x  =  -eB  u/m.c,  (lb) 

where  we  have  chosen  a  Cartesian  coordinate  system  (x,  y,  z)  with  z  along  the 
direction  of  the  magnetic  field  Bq  and  x  along  the  direction  of  wave  proDaaation, 
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(u,  v)  are  the  components  of  the  ion  velocity  in  the  (x,  y)  directions,  respec¬ 
tively,  m.  is  the  ion  mass,  t  the  elapsed  time,  and  c  the  speed  of  liqnt. 

To  mmol  Pt-P  tho  at  ion  .  WO  nroH  in1'  r  '■"I- ' m  i*  t  •>  *  '.>p  . 

Tni/?t  +  S(n.u)/Sx  =  0.  (2) 

We  are  interested  in  travelinq  wave  solutions  of  the  form:  F(x-u  t), 
where  u,  is  the  phase  velocity.  For  such  solutions,  (1),  (2)  and  (i),  (iii) 
can  be  combined  to  yield  a  coupled  set  of  nonlinear  ordinary  differential 
equations : 

^7  =  M2N3  v/(N2  -  M2),  (3a) 

■j -v  =  N-)  ,  (3b) 

where  N  =  n/n  .  M  =  u./c  ,  v  =  v/u.,  f.  =  ft.(x/u.  -t),  U.  =  eB  /m.c  is  the  ion 

0  <?  S  <p  '  }  <p  ^  01 

1 /2 

cyclotron  frequency,  and  c^  =  (kTg/m. )  is  the  ion  sound  velocity.  Without 
loss  of  generality,  we  take  u^  >  o. 

C.  ANALYTIC  SOLUTIONS 

For  supersonic  propagations,  we  have  M  -  I .  The  solutions  of  (3)  can  be 
discussed  conveniently  in  terms  of  the  integral  curves  in  the  v-N  phase  plane 
(Fiq.  la).  Dividing  (3b)  by  (3a)  gives: 

dv/dN  (N-l  )(N2-M2)/M2N3v  'o'- 

Ti  .  ><■>  two  finite  singular  points  for  N  n  1  oca ^rd  a  *  fa'  v  •  N  -  1  . 

era  v  -  o,  N  =  M.  (a)  is  a  center  and  (b)  is  a  saddle  point.  We  note 
i‘  .  p-  ■  iodic  solutions  are  possible  for  o  N  M  The  full  fan>- 
■  •  n*  VPS  m,.y  be  obtained  exp i  ic  ■  r  1  y  '.y  ;:iU-ur  .tuui  (4  .  ,»■ 
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v2/2  +  V  =  K,  (5a) 

with 

V(N)  =  [(1-2N)/2N2  +  (HnN-N)/M2].  (5b) 

The  general  shape  of  V(N)  for  M  >  1  is  depicted  in  Fig.  (lb).  There  is  a  local 
minimum  V,  at  N  =  1  and  a  local  maximum  Vu  at  N  =  M. 

\  n 

Thus,  the  results  of  our  calculations  for  supersonic  propagations  have  the 
simple  analogy  of  a  "particle"  in  a  nonlinear  "potential",  V(N).  The  "kinetic 
energy"  of  the  particle  is  v  /2.  The  "total  energy"  K  parameterizes  the  integral 
curves  in  Fig.  (la).  For  V-|  <  K  <  V^,  the  solution  curves  form  closed  loops. 

Since  dN/d£  >  o  (<  o)  for  v  <  o  (>  o)  when  M  >  N,  the  direction  of  "time"  £, 
circulates  around  the  loops  and  the  corresponding  solutions  for  N(£)  are  there¬ 
fore  periodic.  For  K  >  V^,  the  solutions  are  unbounded.  The  "separatrix  loop", 
which  divides  the  periodic  solutions  from  the  unbounded  solutions,  is  parameterized 
by  K  =  V^.  There  are  no  solutions  for  K  <  .  For  K  =  V-j ,  the  plasma  is 

unperturbed. 

We  now  complete  the  solution  for  M  >  1  by  integrating  (3a)  with  v(M)  prescribed 
by  (5a, b) . 

£  =  /  dN(N2-M2)/M2N3  v(N)  +  C.  (6) 

We  note  that  along  each  branch  (v  <  o  or  >  o)  of  a  closed  loop  with  M  >  N 
and  V-j  <  K  <  VM,  the  integrand  of  (6)  does  not  change  sign.  Thus,  the  integrals 
given  by  (6)  along  the  closed  loops  represent  the  density  profiles  N(£)  for 
periodic  nonlinear  supersonic  propagations.  Since  dN/d£  ~  constant  near  the 
separatrix  loop,  there  are  no  soliton-like  solutions  for  M  >  1. 
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The  electric  field  may  be  expressed  in  terms  of  N. 

E  -  N_1  dN/dC  =  M2N2v(N)/(N2-M2),  (7) 

where  again  v(N)  is  given  by  (5).  Equations  lb;  \  /  constitute  „..e  comp  o.c 

analytic  solution  for  N(£)  and  E(E)  for  supersonic  nonlinear  EIC  waves.  The 

periods  and  waveforms  of  these  nonlinear  waves  depend  on  the  propagating  speed 

(uj  and  the  amplitude  (N  or  K). 

ma  x 

One  class  of  interesting  limiting  solutions  is  obtained  for  large  values 

of  M  and  N  <M.  We  have  plotted  the  profiles  of  N(0  and  E(s)  for  such  a 
max 

case  for  M  =  5  in  Fig.  (2a).  The  electric  field  profile  takes  on  a  distinct 
double-spiked  structure.  For  smaller  values  of  M  and  N  <  M,  the  electric 
field  profiles  vary  gradually  with  £  and  become  sawtooth-like.  For  ,N-1;«1, 
the  solutions  reduce  to  sinusoidal  oscillations  and  the  dependence  on  M  and  K 
disappears.  These  profiles  bear  the  signatures  of  the  linear  and  nonlinear  £  I C 
waves  observed  in  the  suprauroral  region  of  the  maanetosphere  (Mozer  et  ai.,  1 °po ) 
as  was  noted  previously  by  direct  numerical  (Temerin  et  al.,  1979)  and  small  hi/, 
finite  amplitude  calculations  (Chaturveoi,  1976). 

D .  APPROXIMATE  EXPLICIT  SOLUTIONS 

In  fact,  for  sufficiently  large  values  of  K,  (6)  may  be  closely  approximate; 

2 

by  dropping  the  (tnN-N)/M  term  and  the  integral  can  be  carried  out  explicitly. 

For  example,  for  v  -  0, 

£  =  [(2N-1  )/N2  +  2K]1/2  «■  sin^fd-N'1  )/<  i  J 

-  N  1  -  l  1  ^  S  in  V  ( 2KN+  1 }/  (  ;2:  . ) 


^  ■■  mi  1 1  inn  ai. ; 
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For  hypersonic  propagations,  M  »  N  >1,  the  third  term  in  (8)  is  also 

m3  x 

negligible.  In  this  limit, (7)  becomes: 

E-N  [2KN2  +  2N-1]1/2.  (0) 

Thus,  the  profiles  become  independent  of  the  Mach  number  M  explicitly. 

There  remains,  of  course,  the  implicit  dependence  on  the  speed  uf  throuoh  the 

v 

definition  of  £.  The  effective  potential  under  this  approximation  is  characterized 

-1  2 

by  the  competition  of  an  attractive  term  -N  and  a  repulsive  term  1/2N  .  There 
is  no  longer  a  local  maximum  at  finite  N  and  the  saddle  point  has  been  moved  to 
infinity.  We  note  that  in  this  limit,  the  period  of  propagation  is  always  2~ 
units  of  £  and  the  nonlinear  profiles  N(£)  and  E(£)  now  depend  on  the  amplitude 
Nmax  alone.  We  give  three  representative  electric  field  profiles  in  Fig.  (2b) 
for  N  w  =  1.25,  2.5,  and  5.  As  N  „  is  increased,  the  profiles  are  steepened 
exhibiting  sawtooth  and  double-spiked  structures . 
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IV.  EIC  CAVITONS ,  PERIODIC  DENSITY  DIPS,  AND  "PAIRED  SHOCKS"  IN  THE  SUPRAURORAL 

REGION _ _ _ 

A.  INTRODUCTION 

Electric  field  measurements  in  the  suprauroral  region  obtained  from  the  S3-3 
satellite  revealed  the  existence  of  large  paired  electric  fields  ("shocks") 
oriented  predominantly  perpendicular  to  the  terrestrial  magnetic  field  (Mozer 
et  al.,  1980).  Plasma  density-dips  were  sometimes  observed  to  occur  in  conjunc¬ 
tion  with  the  paired  "dc"  electric  fields  which  were  generally  interspersed  with 
coherent  spectra  of  electrostatic  ion  cyclotron  (EIC)  waves. 

In  Sec.  Ill,  we  have  demonstrated  that  double-spiked  and  sawtooth-like  EIC 
waves  are  admissible  solutions  to  the  basic  EIC  equations  for  supersonic  (super¬ 
ion-acoustic)  propagations,  verifying  the  earlier  numerical  calculations  of 
Temerin  et  al .  (1979).  We  now  show  (Chang,  1982)  that  these  same  equations  also 
admit  exact  nonlinear  caviton-like  and  periodic  density-dip  solutions  propagating 
at  subsonic  (sub-ion-acoustic)  speeds.  The  corresponding  electric  field  profiles 
exhibit  paired  structures  with  very  large  (infinite)  peak  intensities.  These 
waves  may  be  interspersed  with  each  other  and  modified  coherently  by  low  amplitude 
EIC  turbulence. 

B.  INTEGRAL  CURVES  FOR  SUBSONIC  PROPAGATIONS 

The  EIC  equations  (equations  3  and  4  of  Sec.  Ill)  are  still  valid  for  sub¬ 
sonic  (M  <  1)  propagations.  Thus,  the  integral  curves  have  the  same  general 
shapes  as  those  for  supersonic  propagations.  In  the  v-N  phase  plane  for  N  >  o 
(Fig.  1),  we  again  have  two  singular  points:  (a)  v  =  o,  N  =  1  (saddle  point), 
and  (b)  v  =  o,  N  =  M  (center).  The  integral  curves  can  form  closed  loops  for 
o  <  N  <  1.  They  are  parameterized  by  V^  <  K  <  V^  and  represent  periodic 
density-dip  propagating  waveforms. 
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The  separatrix  loop  which  passes  through  the  saddle  point  is  given  by 
K  =  .  It  is  easily  seen  from  equations  (5a,  b)  of  Sec.  Ill  that  the  expression 

V(N)  -  V,  -  o  has  a  double  root  at  the  saddle  point-  for  m  <•  1  TH.ir.  the  solu¬ 

tion  correspondi ng  to  the  separatrix  loop  has  an  infinite  period  of  propagation 
and  therefore  has  a  caviton-like  structure. 

C.  DENSITY  AND  ELECTRIC  FIELD  PROFILES 

We  obtain  the  density  profile  N(£)  by  integrating  equations  (3a)  and  (5) 
of  Sec.  III.  We  note  that  dN/d£  is  infinite  at  N  =  M.  Thus,  £  takes  on  extremum 
values  along  the  line  N  =  M.  This  means  that  the  path  of  integration  along  each 
loop  must  jump  from  one  branch  (v  <  o)  to  another  (v  >  o)  at  N  =  M.  Referring 

to  Fig.  1,  two  of  such  admissible  path  of  integration  are  indicated  by: 

(s)  A+3+C'B~C+A~  (for  the  caviton-profile) ,  and  (b)  P+Q+R~Q~R+P~  (for  a  typical 
periodic  profile).  Therefore, 

£  =  /  dN|(N2-M2)/M2N3  v(N) |  +  C  (1) 

along  the  path  of  integration  indicated. 

The  electric  field  E  can  be  calculated  from  equation  (7)  of  Sec.  Ill  We 
note  that  it  is  infinite  at  N  =  M.  At  N  =  M,  the  value  of  v  change*  sign  .v. 
dictated  by  the  condition  that  £  must  increase  (or  decrease)  continuously.  Thus, 
a  vorticity  layer  exists  within  a  thin  sheet  £(M)  -  A£  <  £  <  £(M)  +  A£,  where 
the  kinetic  and  finite  gyroradius  effects  must  be  taken  into  consideration. 

This  kinetic  layer  will  reduce  the  peak  electric  fielc  at  N  -  M  from  infinite 
t'  (but.  still  large)  values.  Previous  calculations  (Yu,  l07?)  worked 

’•'I  ri  £-  density  profile  alone  and  therefore  were  not  aware  of  the  jump  of  the 


7'-' 


i 1 component  v  at  N  =  M. 
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Exact  density  N(0  and  electric  field  E(£)  profiles  for  a  typical  caviton 
are  shown  in  Fig.  (2).  The  electric  field  profile  resembles  those  detected  in 
the  suprauroral  region  of  the  polar  magnetosphere  as  "paired  dc  electrostatic 
shocks"  (Mozer  et  al.,  1980).  The  dominant  feature  of  the  caviton  structure 
generally  spans  over  a  characteristic  distance  of  several  units  of  u^/i^. . 

Periodic  density-dip  profiles  have  similar  forms  except  that  they  have 
finite  periods  of  propagation. 

The  fact  that  the  solution  path  jumps  at  N  =  M  indicates  that  an  actual 
solution  could  be  several  admissible  solutions  interspersed  among  themselves. 

When  plasma  wave  turbulence  (which  is  observed  to  co-exist  with  the  paired  shocks) 
and  other  kinetic  instabilities  are  present,  these  dissipative  effects  can  switch 
the  path  of  integration  to  smaller  and  smaller  loops  eventually  creating  a  density 
cavity  behind  the  wave,  thereby  forming  a  bonafide  shock  structure. 
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of  increasing 


» 


V-l 


V.  CHARGED  BEAM-PLASMA  WAVE  INTERACTIONS 
A .  INTRODUCTION 

We  are  concerned  with  the  generation  of  plasma  wave  radiations  by  a  beam 
of  charged  particles  (electrons  or  ions)  artificially  injected  in  the  ionosphere. 
Such  beam  induced  plasma  wave  experiments  have  been  carried  out  in  space  by 
Cartwright  and  Kellogg  (1974),  the  French-Russ ian  ARAKS  team  (Delahaye  et  al . , 
1978;  Cambou  et  al.,  1978),  and  others  (e.g.,  Haerendel  and  Sagdeev,  1981;  and 
Moore  et  al . ,  1982).  Detailed  discussions  of  these  and  other  related  experiments 
can  be  found  in  a  review  article  by  Winkler  (1Q80). 

Alekhin  et  al .  (1971)  and  Gendrin  (1974)  have  shown  that  for  an  electron 
beam  (typically  10-40  keV,  1-100  ma)  injected  along  the  magnetic  field  lines, 
the  space  charge  effects  and  magnetic  field  strength  limit  the  beam  radius  and 
density  during  the  initial  expansion  of  the  beam  pulse.  This  idea  of  balancing 
the  Lorentz  force  with  the  electric  force  in  vacuum  was  suggested  years  ago  by 
Brillouin  (1945).  The  resulting  differential  equation  admits  both  a  solution 
with  a  limiting  "Brillouin  radius,  r^"  described  above  and  an  oscillatory  solution 
around  rg. 

For  injections  oblique  to  the  magnetic  field,  Jones  and  Kellooq  (1073)  took 
a  simple  modei  of  a  filled  cyclinder.  The  radius  of  the  cylinder  was  taken  ar 
the  Larmor  radius  r^  =  ux/fiQ,  where  ux  is  the  perpendicul ar  velocity  of  the 
beam,  and  the  electron  cyclotron  frequency.  The  lenqth  of  the  cylinder  was 
estimated  from  the  combined  diffusion  effects  produced  by  the  enerqy  soread  ot 
the  electrons,  the  angular  dispersion  of  the  beam,  and  the  finite  duration  of 
the  pulse  (-1-20  ms)  using  simple  adiabatic  invariant  arnuments.  With  a  warm 
ionospheric  plasma  background,  they  showed  that  a  variety  of  electrostatic  and 
electromagnetic  plasma  waves  can  be  excited.  They  discussed,  in  particular. 


the  electrostatic  plasma  mode  of  a  band  of  radiations  between  the  electron 
plasma  frequency  and  the  upper  hybrid  frequency  Juh ,  the  electrostatic 
Bernstein  modes  around  n'i  .  and  the  elerf  ror'anno*-  1  r  whistler  mr>apc;  below 
Assuming  tne  distributions  are  spatially  uniform,  tney  calculated  tne  growtr 
rates  and  the  parallel  and  perpendicular  components  of  the  group  velocities  of 
the  plasma  and  whistler  modes.  Using  these  results  and  tne  assumed  "filled 
cylinder"  beam  geometry,  they  made  rough  "order  of  magnitude"  estimates  of  the 
waves  amplitudes  and  radiative  powers  of  the  calculated  unstable  modes. 

Lavergnat  and  Pell  at  (1979),  on  the  other  hand,  used  both  a  hollow  cylinder 
(radius  r,  ,  tnickness  rD),  and  a  more  realistic  helix  model  (helix  radius  r.  . 
beam  radius  rg)  to  calculate  the  plasma  waves  spontaneously  produced  (e.g.,  via 
Cerenkov  and  cyclotron  radiations)  by  the  oDlique  electron  pulses.  It  was  shown 
that  plasma  modes,  Bernstein  modes,  and  whistler  modes  can  again  be  excited. 

Both  electrostatic  and  electromagnetic  plasma  wave  radiations  generated  by 

artificially  injected  electron  beams  in  the  ionospnere  have  been  observed  experl 

mentally.  These  include  the  plasma  waves  discussed  above  as  well  as  the  electro 

static  "lower  hynno  waves  and  the  predomi  nantly  electromagnetic  VLB  wave:,.  Tf-' 

frequency  range  of  these  waves  are  characterized  bv  .1 .  <«.  <•:  >  and  their 

i  e 

wave-vet  tor?  k  generally  make  large  oblique  angles  with  respect  in  .he  amoi 
ignetic  field.  Ihese  waves  couple  the  dynamics  of  the  electron  beam  and  the 
amniort  ions.  In  fact,  they  can  effectively  energize  the  ions  in  the  transverse 
d  .•  •.*'  .•  witn  respect  to  the  ambient  magnetic  f  i  el  d  ires.  c.xpf  riaer/a1  •  • 

•  ‘  such  ivpe  cf  transverse  ion  acceleration  have  been  reported  c.v  Arnoi.iv 

-  r  :  i  c.pg  , . 

Inn  beam  injections  in  the  ionosphere  are  very  different  from  those  of" 

•  '  "earns  .  Firstly,  we  expect  the  ion  beams  to  srr-\ad  i  n :  n  > "  '  v  '  '■.<> '  i 
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larger  dimensions  (because  of  the  larger  Brillouin  and  Larmor  radii)  even  with¬ 
out  wave-particle  interactions.  This,  in  turn,  allows  the  unstable  modes  to 
interact  at  longer  time  intervals  and  therefore  possibly  attain  larger  wave 
amplifications.  Secondly,  new  unstable  modes  are  expected  from  ion  beam  injec¬ 
tions.  These  include  the  electrostatic  ion  cyclotron  waves,  non-resonantly  excited 
ion  cyclotron  harmonics,  ion  acoustic  waves,  and  the  broad-band  electrostatic  and 
electromagnetic  VLF  waves.  These  waves  can  interact  with  the  ion  beam  causina  it 
to  spread  anomalously.  Furthermore,  the  low  frequency  electrostatic  waves  can 
parametrically  interact  with  the  high  frequency  waves  in  the  magnetosphere  and 
convert  them  to  el ectromagnetic  radiations.  Thirdly,  ion  beams  are  known  to 
spread  anomalously  in  an  ambient  plasma  containing  pre-existing  fluid  turbulence 
(e.g.,  those  due  to  Kel vin-Helmholtz  instability)  when  electrostatic  ion  cyclotron 
waves  are  also  present.  Finally,  as  it  will  be  shown  presently,  lower  hybrid 
waves  excited  by  the  ion  beam  can  accelerate  the  electrons  in  the  field-aligned 
direction.  The  energized  electrons  can,  then,  excite  other  secondary  plasma 
waves . 

A  series  of  interesting  laboratory  experiments  designed  to  simulate  actual 
electron  beam  injections  in  space  environments  have  been  conducted  by  Bernstein 
and  his  co-workers  (Berstein  et  al.,  1975,  1978,  and  1979;  Denig,  1982).  The 
electron  gun  was  generally  operated  dc,  although  sometimes  a  pulsed  gun  was  also 
used.  For  a  1  keV  beam  with  I  <  10  ma  at  10"^  torr  in  a  magnetic  field  B  =  1.5  G, 
the  beam  assumed  a  multi-node  configuration  with  a  node-length  where 

u„  is  the  parallel  beam  velocity.  RF  measurements  indicate  fluctuations  near 
~  and  the  lower  hybrid  frequency. 

When  the  current  is  raised  (at  a  fixed  beam  energy)  beyond  some  critical 
value,  the  RF  fluctuations  form  (i)  a  band  of  emissions  near  the  plasma  frequency 
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and/or  upper  hybrid  resonances,  (ii)  a  board  band  of  whistler  mode  emissions 
with  uj  <  3  ,  and  (iii)  an  intense  broad  band  of  VLF  emissions.  A  column  of 
plasma  of  enhanced  density  is  observed  and  there  is  evidence  that  the  primary 
electron  beam  is  dispersed  both  in  energy  and  Di ten  angle.  Ihe  peax  luminosity 
of  the  discharge  becomes  about  ten  times  that  of  the  beam  alone  and  the  radius 
is  increased  nearly  three-fold. 

Such  type  of  beam-plasma-discharges  (BPD)  can  be  important  in  actual  space 
experiments.  The  increased  plasma  density  can  produce  the  observed  neutral i cation 
of  the  space  vehicles  by  providing  the  required  return  current.  The  process  is 
accompanied  by  strong  wave-particle  heating,  severe  energy  diffusion,  pitch  angle 
diffusion,  and  radial  expansion  of  the  beam  profile,  placing  practical  limitations 
on  actual  beam  injections  in  space. 

It  nas  been  recognized  that  plasma  instabilities  and  me  ensuino  wave- 
oarticle  interactions  are  required  to  accelerate  the  cnarueu  nartie'es  in  order- 
re  Produce  the  observed  anomalous  ionizafon  in  BPd  experiments.  Easier  rneor*-vca 
werx  ( Papadopoul os ,  1932)  is  based  on  the  high  frequency  electron  modes  driven 
unstable  by  the  electron  beam.  It  is  our  contention  : C-nnc  and  asperse,  ; 3  v ' 
that,  the  low  frequency  branch  of  the  beam-excited  plasma  wave-  can  be  very  'Vncrrant 
in  determining  the  coupled  electron  and  icn  dynamics  of  3PD.  We  suggest  tna’  ir. 
me  presence  of  an  external  magnetic  field,  lower  hybrid  waves  mav  play  an  im.'crmn: 
role  in  determining  the  critical  current  for  the  onnet  of  BPD.  In  addition. 

•y  r  i  -k---  whist'ers  (or  electromagnetic  lower  hybrid  waves'  are  i'  ?*■  ire  -,. 

^  'osta;1';  lower  hybrid  waves  carry  most  of  the  energy  and  can  r-esonjit  , 

.,e  ^  particles  tc  the  ionization  ■  >-eshol  1  of  the  r.rpV-  •  -  -  , 


r q  fro  enhanced  ionization  detected  in  BPP  experiments. 
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B.  LOWER  HYBRID  WAVES  AND  HYBRID  WHISTLERS  PRODUCED  BY  AN  ELECTRON  BEAM 

It  is  known  that  in  a  magnetoplasma,  low  frequency  ( oj  «  L2e)  whistler 

modes  can  be  excited  by  an  energetic  electron  beam.  The  modes  that  propaqate 

at  very  large  oblique  angles  with  respect  to  the  magnetic  field  are  excited  by 

the  electron  beam  primarily  through  Landau  resonance.  These  modes  (with  k„  <<  kA  ) 

can  transfer  energy  to  the  ions  in  the  transverse  direction  via  ion  Landau-  and 

gyro- resonances .  To  simplify  the  calculation,  let  us  search  for  modes  with 

?i.  «  ^  «  ii  and  o-  >>  kx  ^  >>  p  .  [w  <  \1.  can  also  exist.  We  will  return  to 
i  e  i  e  i 

this  point  later.]  Thus,  we  work  under  the  approximation  that  the  electrons  are 
strongly  magnetized  and  the  ions  nearly  unmagnetized. 

The  dispersion  relation  for  >>  is  then  approximately  (Marsch  and 
Chang,  1982): 


D(u)  =  a 2^  "  a]  ^  +  a0  =  o. 


0) 


where 


ao  =  ('VV2  ma2Cma2  '  (1+k„2/k2)3(  1  -Ze) 


(2a) 


a!  =  (k„2/k2)S2e2[l-Ze-(a)i/coe)2Zi] 

+  (k,2/k2)|nen.|(i-zi) 

+  MA2|neQi|[(l+k„2/k2)Ze-2]  (2b) 

a2  =  l-(k„2/k2)Ze  -  (oii/u>e)2Z1 ,  (2c) 
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with  Mft2  =  (oo/kVA)2, 

VA2  =  B2/E4)rnj.mj, 

Z  =  E(n./n  )[l+(oj/kllv.„)tWlw  /kl(v.(1jj, 
e  j  j  e  j  j 

Zi  =  Z(to^/u)i ) 2 [  1  +( lo/ k v^( ot )  )2W(co'  /  kv^(cx) )], 

a 

2 

W(z)=  (2*)'1/2/°dx  x  e'x/2/(x-z),  (3) 

to  =  to-k„u  j  , 

2 

V-  H  ■“  k_T  ,/m., 

J  ts  J  J 

2(a)  =  Vi2„  cos2oi  +  sin2a, 

2  "  2  2  2  2 
cos  a  =  ki.^/k  ,  sin  a  -  kj.  /k  . 
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(yMr)a  -  Im  D(w) , 


where 


Im  D(<2)  = 


5j2(<22-|fi  !). 


.[(k„/k)2  Im  Zg  +  (ui/ufe)t  Im 

-  r((22-l)(l+k„2/k2)f1A2  +  MA4]Im  Z£ 


We  can  conclude  from  (5)  and  (6)  that: 


(i)  the  electromagnetic  corrections  are  small  on  the  resonance  cone, 

(ii)  ion  damping  (heating)  is  important  for  near  perpendicular 

2  2 

propagations  (k„  /kx  -  n^/uK ) , 

(iii)  the  growth  rate  y  is  large  on  the  resonance  cone. 

C.  THE  ROLE  OF  LH  WAVES  IN  BEAM- PLASMA- DISCHARGES 

From  the  discussion  of  Sec.  VB  ,  we  expect  that  both  the  electrostatic  LH 
waves  (electrostatic  whistlers)  and  the  electromagnetic  LH  waves  (i.e.,  "Hybrid" 
or  very  obliquely  propagating  whislters)  will  be  excited  by  an  energetic  electron 
beam,  with  the  former  carrying  most  of  the  energy. 

The  primary  electron  beam  is  expected  to  ionize  initially  a  small  fraction 
of  the  ambient  medium  and  at  the  same  time  to  generate  the  LH  waves.  The  electro¬ 
static  LH  waves  transfer  the  energy  of  the  electrons  to  the  ions  [in  fact,  energiz¬ 
ing  a  fraction  of  the  ions  above  the  thermal  energy].  When  some  of  the  ions  are 
energized  above  the  ionization  threshold,  more  of  the  ambient  medium  is  ionized. 
This  process  (indicated  by  the  solid  lines  of  the  schematic  loop  diagram  in  Fig.  1) 


? 
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continues  until  that  part  of  the  lower  hybrid  wave  energy  responsible  for 
energizing  the  ions  just  above  the  ionization  potential  is  saturated.  If  the 
electron  beam  is  maintained,  the  inni^  +  i^n  viH  ™^h  '  f 

state.  The  total  amount  of  ionization  in  this  state  is  expected  to  be  much 
higher  than  that  can  be  produced  without  the  LH  wave-particle  interaction. 

This  process  of  anomalous  ionization  is  at  least  partially  responsible  for  the 
phenomenon  of  BPD. 

Recent  BPD  experiments  in  the  large  vacuum  chamber  at  the  Johnson  Space 
Center  carried  out  with  the  Plasma  Diagnostic  Package  (PDP)  designed  for  the 
space  shuttle  (Denig,  1982)  indicated  the  existence  of  a  trapped  ion  population 
(typically  with  an  average  energy  of  50-60  eV)  upon  the  onset  of  BPD.  There  was 
also  the  indication  of  the  appearance  of  a  radial  dc  electric  field  with  a  potential 
well  located  near  the  beam  gyroradius.  These  results  are  tentative  and  require 
further  experimental  verification  based  upon  more  refined  diagnostic  techniques 
and  instrumentation.  Furthermore,  the  theoretical  justification  of  this  kind 
of  ion  trapping  does  not.  appear  to  be  immediately  apparent.  If  the  trapping  does 
set  in,  however,  it  can  be  argued  that  the  positive  pei  pendic.ular  gradient  cf 
the  resulting  ion  distribution  function  can  excite  the  ion  cyclotron  harmonic 
waves  (Bohmer,  1976;  Kintner,  1980)  as  well  as  secondary  lower  hybrid  waves. 

Both  of  these  waves  can  effectively  energize  the  electrons  in  the  field-aligned 
direction.  These  energized  electrons  can,  in  turn,  ionize  the  ambient  medium 
uni  tnu  newly  produced  low  energy  electrons  can  again  be  accelerated  by  the 
seionJ.  ry  LH  waves  and  the  ion  cyclotron  harmonics.  Beyond  a  certain  critical 
■  .  .  this  effect  can  efficiently  multiply  to  produce  a  s  ion5 ;  «<  add ;  ‘  v  m  1 

f>  tI ion  of  the  enhanced  ionization.  Such  a  process  may  be  quite  temporal  ’nd 
v  ■  I  * 1  this  may  provide  an  explanation  of  thrt  c  b-^rv'-d  '  i.'  hf  at  i  v 
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of  the  total  intensity  of  ionization.  We  have  indicated  this  "probable"  process 

as  broken  lines  in  the  loop  diagram  of  Fig.  1. 

We  note  that  a  similar  process  also  occurs  for  a  high  current,  energetic 
ion  beam  (see  Fig.  2).  Such  processes  may  have  been  observed  by  the  PORCUPINE 
Xe+-gun  experiment  (Kintner,  private  communication),  and  Ar+-nun  sounding  rocket 
experiments  (Moore  et  al.,  1982). 
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VJ.  THE  UBIQUITOUS  LOWER  HYBRID  MODES  IN  SPACE  PLASMAS 

A.  INTRODUCTION 

We  have  seen  that  the  electrostatic  lower  hybrid  waves  (electrostatic 
whistlers)  and  electromagnetic  lower  hybrid  waves  (hybrid  or  very  obliquely 
propagating  whistlers)  populate  the  topside  ionosphere,  the  suprauroral  region 
and  are  detected  in  naturally  occurring  or  artificially  created  beam-plasma 
processes.  Our  investigations  have  shown  that  lower  hybrid  waves  are  omnipresent 
in  the  ionospheres  and  magnetospheres  of  the  Earth  and  other  planets,  ahead  of 
the  planetary  bow  shocks,  across  interplanetary  shocks,  and  in  the  solar  wind. 

Space  does  not  allow  a  full  discussion  of  these  and  other  phenomena  involving 
lower  hybrid  wave-particle  interactions  in  space  plasmas.  As  an  example,  we 
will  demonstrate  that  the  frequently  observed  broad-band  low  frequency  electro¬ 
static  noise  in  the  solar  wind  in  fact  has  a  dominant  lower  hybrid  component. 

Other  similar  examples  are  abound  and  have  been  detected  in  plasmas  of  various 
terrestrial  and  extra-terrestrial  environments. 

B .  LOWER  HYBRID  WAVES  IN  THE  SOLAR  WIND  (Marsch  and  Chano,  1982) 

Plasma  wave  measurements  in  the  solar  wind  (Gurnett  and  Frank,  1978;  Gurnett, 
1981;  Scarf  et  al.,  1981)  often  exhibit  broad-band  low  frequency  electrostatic 
turbulent  noises  accompanied  by  bursts  of  low  frequency,  obliquely  propagating, 
hybrid  whistler  mode  fluctuations  (Beinroth  and  Neubauer,  1981;  Coroniti  et  al . , 
1982).  These  electrostatic  noises  are  usually  identified  as  ion  acoustic  waves 
(Gurnett  and  Anderson,  1977;  Gurnett  et  al.,  1979).  However,  because  of  ion 
damping,  ion  acoustic  waves  can  be  excited  by  the  solar  wind  particle  distribu¬ 
tions  only  when  the  electron  temperature  Tg  is  much  larger  than  the  ion  temperature 
T. .  Most  of  the  time  this  condition  is  not  satisfied  when  the  low  frequency 
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electrostatic  fluctuations  are  detected  ( Dum  et  al.,  1981).  Thus,  an  alterna¬ 
tive  explanation  of  these  frequently  observed  e lectros tatic  modes  is  needed. 

We  propose  here  that  lower  hybrid  wave0  rap  hr  l-,'-,  r  r  +  -i -- 

tron  distributions  that  represent  the  heat  flux  in  the  solar  wind.  These  modes, 
Doppler-shifted,  are  expected  to  have  the  observed  frequency  range  of  the  LF 
electrostatic  noises. 

C .  DISPERSION  RELATION 

The  solar  wind  electron  heat  flux  distribution  can  be  modeled  by  a  thermal 
core  and  a  suprathermal  halo  component  (Feldman  et  al.,  19751.  It  is  known  'hat 
such  an  electron  distribution  can  excite  electromacinetic  wnistler  modes  (Manns. 

1976).  Neglecting  electron  damping  of  the  core  distribution  for  tiie  moment,  the 

2  2 

modes  that  cluster  around  the  lower  hybrid  resonance  cone  for  k,,6  -  kx  and  a 

cold  plasma  obey  the  dispersion  relation 

.o2  =  ill  .2[l+(k„i»>  /kui  .)^j/[l+(m  / T!  ) ^  1 ,  (1' 

pi  L  '  ’  pe  pi'  J  1  '  pe  e' 

>2  2  2 

whore  .  =  4irn  e.  /tit.,  and  m _  =  Tnn  e  /m  are  the  squares  of  the  ion  and 

pi  ill  pe  e  e 

electron  plasma  frequencies,  11  is  the  electron  gyrofrequency ,  and  (k„,kx)  are 
the  parallel  and  perpendicular  component  of  the  wave  vector,  k.  For  . 

this  expression  reduces  to  that  discussed  by  Coppi  et  al.  ( 1 0 76 )  and  Papadopoulos 
and  Palmadesso  (1976)  in  the  electrostatic  limit  and  the  magnitude  of  the  frequency 
is  of  the  order  of  the  ion  plasma  frequency.  In  the  solar  wind,  on  the  other  nan.,!, 

.  ‘  e" ,  and  (1)  gives  .>  -  jll  Sh|  ,t,he  classical  lower  hybrid  >-esonance 
t  i ’-o  -rv.  Thus,  the  magnitude  of  the  frequency  re  nor  of  the  lowe*-  hybrid  "•odes 

> '■  •  solar  wind  can  be  much  less  than  that  of  the  ion  plasma  frequency  with 
■  ■  ■  ,  where  i>.  is  the  ion  gyrofreouenny . 
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In  this  discussion,  we  consider  purely  electrostatic  modes  in  an  electron- 
ion  plasma,  where  the  ions  can  be  treated  as  unmagnetized.  The  full  dispersion 
relation  may  be  written  as 

k2  =  -  z  tj  Ej  (k,u>),  (2) 

2 

where  tj  =  (kj/k.)  is  the  squared  ratio  of  Debye  wave  vector  of  the  jth  species 
2  2 

(k.  =  4tt  e,  n  ./T  .)  to  the  ion  Debye  vector  k.  and  k  =  k/k . .  In  the  unmagnetized 

J  J  J  J  1  1 

case,  the  normalized  dielectric  constant  for  a  bi-Maxwel 1 ian  with  a  mean  squared 
2  2  2 

thermal  speed  v.  =  (Vj„  +  2Vjj.  )/3  is  defined  by 

Cj(k.oj)  =  (V/Vj(a)]2  W[oj'/kVj(a)],  (3) 

with  the  Doppler-shifted  frequency  co ’  =  to  -  k„u,.  and  the  effective  thermal  velocity 

J 

2  2  2  2  2 

v-  (a)  =  v cos  a  +  v.x  sin  a.  The  function  W(z)  is  related  to  the  olasma 

J  J  »J 

dispersion  function  (Ichimaru,  1973).  The  angle  a  determines  the  wave  propagation 
direction  (k„  =  k  cos  a  and  kx  =  ksina)  with  respect  to  the  maaneti  c  field.  We 
allow  for  streaming  along  the  field  relative  to  the  center  of  mass  frame  with 
velocity  u^  for  species  j.  The  wave  phase  velocity  refers  to  the  same  frame. 

The  magnetized  dielectric  constant  can  be  cast  into  the  form 

£j(k-'“>  '  (YV>2  Ujo  "Jl  Ejn>-  ('" 

The  zeroth  order  contribution  e.  is  the  only  term  that  survives  in  the  stronoly 

jo  J 

magneti zed  limit. 


■;jo  '  1  j  'o(bj>  “  (zo>  -']• 


and  the  gyroresonances  are  described  by 
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in  =  {zo^W(2n^1/zn  *  <H<z-n>-]  )/z-n^ 


-b. 


!  1  ]  T  W  ( z  )  +  W(z  n)]}  e  J  I^b.) 


(4b) 


As  usual,  In(b-)  is  the  modified  Bessel  function  of  order  n  with  argument 


bj  •  ( kj.V'V2  •  i?  (?j‘xj>2- 


(4c) 


For  convenience,  we  nave  introduced  the  normalized  perpendicular  phase  velocity 

=  w/kxv.j.  and  the  frequency  <o  =  cj/Up.  normalized  to  the  ion  plasma  frequency 

The  ration  x.  =  .  determines  the  deqree  of  magnetization  witn  .  =  e.B/m.c. 

J  J  pi  '  J  i  ‘ 

The  normalized  resonant  speed  is  z n  =  ^.,,(1  -  nx^/T)-  u./v.„  whereby  =  ../k„v 

r  Kj 

is  the  parallel  phase  velocity  in  units  of  tne  thermal  speed  'or  tne  1  specie^ 

These  quantities  are  useful  in  localizing  tre  wave-part  ide  resonances  in  veioci 

space.  For  given  parameters  of  tne  particle  d’Stribut ;ons ,  \2)  depends  ;n‘<ca’. 

on  t.  =  (T./T.)  •  (n./n.)  •  (e./e)  and  x-  which  is  determined  by  the  density 
J  1  .1  J  i  J  1  J 

and  the  magnetic  field  strenqth.  Finally,  we  quote  the  normalized  growth  ra~e. 
whicn  for  small  growth  is  qiven  by 


■•'/..)  ■  -  ::  t  Itn  ,  t,  --3Re>:  ./3u.  •: 

j  J  i  1 

0 .  SOME  NUMER yC AL  RESULTS 

We  briefly  discuss  tne  various  possible  low  frequency  modes  derived  from  [2 
>'  ■ -v  •  gy  sources  for  these  modes  are  also  proposed  in  terms  of  :  yOical  solar 
ic  ’  1  :  i c i e  distributions.  The  normalized  dielectric  constants  ,:n  ;4)  .  'n 

*  ■’d  tunc  *  ;  r '  r  3  of  their  arguments  (for  examn  -  W  1  -  • 1  ■:  i  • 

•  n»r.’f  erom  0,02  to  2.0).  For  t  =  1  or  T  =  T.  ,  both  species  yield  comcarab 
'  1  r  ''  i!  inns  t  o  t  ? )  l  f  i  ^  x  =  g  j  „ .  This  ccnrilt  a.n  imol  >es  :k  ~/k  ■■  r  /n  .  . 
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If  : . .  -N1  and  !o  „ '  •  -1 ,  one  has  the  classical  undamped  LH  waves.  Otherwise, 

1  6 

the  dispersion  characteristics  are  modified  by  thermal  effect  and  Landau  dampin'-:. 
Considering  t  "C  !  or  T  one  finds  that  the  electron  contribution  to  (2)  can 

be  kept  small  even  for  '-i  „|<<1  and  ^ j.  - 1  .  Then  one  obtains  the  usual  ion  sound 
mode  of  ion  plasma  oscillations,  Debye-shielded  by  the  background  electrons.  For 
Maxwellian  distributions,  this  mode  can  exist  only  for  t  <  0.285.  Realistic 
distributions  have  been  discussed  extensively  by  Dum  et  al .  (1081). 

Many  plasmas  are  characterized  by  moderate  ratios  of  tg>l.  Besides  the 
classical  lower  hybrid  waves,  another  possible  situation  is  for  !<i  and 
d.jx-'l-  The  last  condition  is  necessary  in  order  to  keep  the  ion  contribution 
low  in  the  total  dielectric  constant.  If  is  large  enough,  the  term  ■  .  in 
(2)  can  be  neglected  and  a  pure  electron  mode  emeroes. 

Retaining  only  terms  to  lowest  order  in  bp  qives  the  following  expression 
for  the  normalized  wave  frequency  in  a  bi-Maxwel 1 ian  electron-ion  plasma  without 
drifts  for  oblique  propagations 

Cl2  =  [~zi2W( zi )  -  (k„/k)2(c:pe/(Jupi.)2  Oe„2WUe„)]/ 

[l+(kx/k)2((ape/ue)2(l-(Tex/Te„)W(^„))],  (6) 

2  2  M2 

where  z^  =  u>/kv.(a)  =  +  $.4.  )  In  the  cold  plasma  limit 

2 

-z  W(z)>l  and  (6)  reduces  to  the  expression  discussed  in  the  introduction  for 
kx  >>  k„.  As  long  as  b  <1,  the  above  result  is  fairly  general  and  can  include 
the  thermal  effects.  For  given  phase  velocities,  the  frequency  is  readily 
obtained  for  (6),  that  yields  a  broad  band  spectrum  in  frequency  which  depends 
on  u)pe/ue,  and  the  wave  propagation  direction.  Solutions  may  also  exist  for 
moderate  values  of  k„/kx<l. 


The  following  numerical  results,  based  on  the  full  dispersion  equation  ,'2;, 

illustrate  the  above  qualitative  arguments.  For  the  sake  of  simplicity,  we  assume 

.3 

Maxwellian  distributions  with  no  =  n?  =  10  cm  ,  v^  =  45  km/'s  and  v^  =  1,250  . 

tnese  are  typical  densities  and  tnermai  sueeas  tci  a  sumr  '»uio  *.  •  .i  vt  iot  i  . ; 

of  about  500  km/s  (Feldman  et  al . ,  1975).  The  temperature  ratio  is  then  "0/'e  =  2.38 

The  ratio  x.  =  ;2./w  .  varies  to  a  certain  extent  in  the  solar  wind  depending  on  tne 
l  l  pi 

macroscopic  stream  structure.  For  x  =  5*10’^,  ik„/kx'  =  1/25,  and  x/k„  =  3,750  km/s 

r 

one  obtains  =  0.054  and  y/u>  =  -0.079.  The  proton  and  electron  damping  are  roughly 
equal  and  the  value  of  w  is  found  to  be  very  close  to  the  value  given  in  (6).  The 
electrons  are  strongly  magnetized  (i.e.,  Laudau  damping  dominates  over  cyclotron 
damping)  with  bp  =  0.24  and  xg  =  0.918.  Actually,  for  the  same  particle  parameters 
used  above,  a  solution  exists  for  various  values  of  xg.  For  example,  for  x  =  13.4 

_  O 

or  x  =  10"1",  the  same  value  of  b  results  if  <L  is  increased  to  0.975.  It  is  worth 
P  ^ 

emphasizing  that,  in  contrast  to  ion  sound,  the  lower  hybrid  waves  do  not  crucially 
depend  on  the  temperature  ratio  Te/Tp. 

In  Table  1,  we  have  listed  some  typical  fast  solar  wind  parameters  with 
T„/T  =  :.  The  electron  distribution  is  represented  bv  two  drifting  bi-Maxwe1 1  ians 
fcore(c)  and  halo(h)].  Using  these  parameters  and  for  w/k„  =  3,500  km/s,  'j/kx  =  464 
km/s,  ip/jpp  -  5-10  a  =  82.4°,  bp  =  167,  we  find  i  =  0.1  and  -/~  =  +0.031.  Leav¬ 
ing  fhe  pnase  velocities  and  other  parameters  unchanged,  an  increase  in  .Tg/ij  to 

18.37  increases  <i>  to  2.0  with  the  same  y.  We  have  chosen  a  value  of  n./n,  =  0.15 

n  p 

-viiirn  is  sligntlv  in  excess  of  what  is  usually  observed.  The  sign  of  •••  is  changed 


■  reducing  the  value  n,  to  0.05  n  . 

^  h  p 

•  mil  1  y ,  we  consider  a  case  with  ii  /,; 

P  PP 


,  a  mii.ihe*'  <•<„.  or. a  i 


at  0.  i  4ij.  For  isotropic  core  and  halo  Maxwellian  distributions  with  vc  =  l,2f)o  ^.m/s . 
r  •  00  km/s,  ur  r  -200  km/s,  uh  =  3,60(1  km/s.  n./n^  -  0.95,  n,  n  -  0.0!,. 


'•094,  ;»i/k„  -  3,500  km/ s  and  iii/k A  ■  328  km/s.  we  find 
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and  x  =  84.6°.  By  slightly  increasing  the  halo  density,  this  wave  can  be  driven 
unstable  by  the  halo  electrons. 

E.  SUMMARY 

We  have  given  numerical  examples  indicating  that  typical  solar  wind  plasmas 

(0.3  -  1  AU)  can  support  electrostatic  lower  hybrid  modes.  Because  to  •  >  -  in 

the  solar  wind,  the  cutoff  frequency  of  these  electrostatic  modes  is  near  the 

1 12 

classical  lower  hybrid  resonance  ~  |  '  .  These  modes  are  probably 

driven  unstable  by  the  anisotropic  halos  (the  heat  flux)  of  the  electron  distribu¬ 
tions.  This  may  have  important  consequences  for  transport  processes. 

Since  these  waves  propagate  nearly  normal  to  the  magnetic  field  lines,  they 
can  energize  the  ions  in  the  solar  wind  in  the  transverse  direction  and  thus  pro¬ 
duce  anisotropic  ion  distributions.  Alternatively,  these  waves  could  be  excited 
by  loss-cone  type  ion  distributions,  that  may  originate  from  interplanetary  shocks. 

The  lower  hybrid  modes  are  generally  accompanied  by  broad  band  electromagnetic 
whistler  modes  in  the  solar  wind  (8einroth  and  Neubauer,  1981;  Coroniti  et  al., 
1982).  These  low  frequency,  very  obliquely  propagating  whistler  modes  can  be 
generated  efficiently  by  the  electron  halos.  We  call  these  waves  the  electro¬ 
magnetic  lower  hybrid  or  hybrid  whistler  modes.  Unlike  the  ion  acoustic  modes, 
the  em  and  es  LH  modes  can  be  supported  for  situations  with  T .  -  Tg  and  may  be 
ubiquitous  in  space  plasmas.  They  have  been  found  in  the  planetary  ionospheres 
and  magnetospheres  (Temerin  et  al.,  1981)  ahead  of  planetary  bow  shocks  (Anderson 
et  al.,  1981)  and  across  interplanetary  shocks  (Coroniti  et  al.,  1Q82). 

The  results  of  the  present  analysis  stronaly  suqaest  that  the  frequently 
observed  low  frequency  electrostatic  noises  and  very  obliquely  propagating  whistlers 
(in  the  solar  wind  and  other  space  plasma  environments)  are  probably  lower  hybrid 
modes  similar  to  those  discussed  in  Sections  II  and  V. 
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VII.  STRONG  PLASMA  TURBULENCE  AND  STOCHASTIC  HEATING 


A.  INTRODUCTION 

In  the  previous  sections,  we  have  discussed  a  number  of  plasma  instabilities 
and  the  associated  wave-particle  interactions  using  the  ideas  of  weak  turbulence 
and  considered  their  possible  relation  to  the  various  observed  ionospheric  and 
magnetospheric  plasma  phenomena.  Most  of  the  observed  plasma  wave  data,  however, 
display  strong  turbulence  signatures.  For  strong  plasma  turbulence,  it  is  well 
known  that  the  correlation  effects  and  response  functions  are  extremely  long- 
ranged,  and  therefore,  cannot  be  adequately  described  by  the  existing  classical 
quasi  linear  or  nonlinear  mode-mode  coupling  plasma  theories.  In  a  strongly 
turbulent  plasma,  charged  particles  can  be  efficiently  accelerated  by  the  lonq- 
ranged  interacting,  stochastic  wave  fields  to  very  large  energies  without  the 
classical  resonant  wave-particle  interactions. 

The  concept  of  large  scale  stochasticity,  strong  plasma  turbulence,  and 
stochastic  heating  and  diffusion  is  still  in  the  intital  developing  stage.  It 
is  the  consensus  that  these  effects  can  be  treated  analytically  using  the  newly 
developed  tool  in  statistical  mechanics--the  renormalization-group.  We  have 
made  a  number  of  recent  contributions  in  this  direction  (Chang,  Ni col  1  and 
Young,  1978;  Vvedensky  and  Chang,  1982;  Vvedensky,  Chang,  and  Ni col  1 ,  1983). 

It  is  shown  that  the  stochastic  description  of  (weak  and  stronq)  plasma  turbulence 
can  be  constructed  in  terms  of  a  path  integral  formalism  similar  to  that  is  employed 
in  quantum  electrodynamics .  Our  starting  point  is  a  generalized  Lanqevin  equation 
based  on  the  Vlasov  formulation  coupled  with  a  random  noise  term.  We  demonstrated 
that  the  "probability  density  functional"  can  be  calculated  using  an  appropriate 
functional  Jacobian  involving  the  noise  fields. 
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To  calculate  the  correlation  and  response  functions,  we  have  constructed 
a  generating  functional  using  the  Legendre  transform  of  the  probability  density 
fucntional .  This  generating  functional  may  be  evaluated  approximately  for  weak 
turbulence  using  tne  quasninear  prescriptions  01  older  truncated  nonlinear 
mode-mode  coupling  interaction  scnemes  provided  the  basic  physical  mechanism  of 
the  plasma  instability  is  known.  Earlier  renormalization  theories  of  Dupree 
(1966)  and  the  direct-interaction  approximation  of  Kraichnan  (1964,  1970)  can  be 
derived  using  this  path  integral  formalism. 


B.  RENORMALIZATION  GROUP  AND  OPERATOR  ALGEBRA  FOR  DYNAMICS 

Based  on  this  path  integral  formulation,  we  have  developed  differential 
renormalization-group  techniques  to  evaluate  the  random  dynamics  of  strongly 
stochastic  nonlinear  systems.  This  renormalization-group  equation  for  dynamics 
can  be  solved  most  conveniently  usinq  an  operator  expansion  technique  which  we 
developed  recently  (Vvedensky  and  Chang,  1982).  Briefly,  the  procedure  is  as 
follows . 

We  start  with  the  one-particle  irreducible  differential  renormalization- 
group  generator  for  dynamics  (Chang,  Nicoll ,  and  Younq,  1979). 


~A 


(d+z)A  +  /- 


dL 

(2  c)< 


Z7!  ‘  V1 


( kw)  [^-(2-n-d)  -k  ■ 


k  Zu  9uj-  6^  i  (  koi) 


+ 


L  9 ( k*, 


[^(2-n-d)-z-k-vk 


-Zu 


3  .  5  A  . 

So;"1  (Si/r^kwT 


+  z  'TVTd  C  2T  tr  ln{A(qu;-q,-o 
\  ^  ) 


r  dp  -  dp  .du1  -dw" 
'(2c)d  ' (2i)  '  " 


-;0_xd  J2=r  f2T  A(q^;p'-')^A  p ' 


"  1 A  {  p "  " ;  -  c  - . ) 

(11 
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where  A  is  the  action  integral,  i^(kw)  and  ^2(kw)  are  the  Fourier  components 
of  the  order  parameter  and  the  conjugate  momentum,  respectively,  with  i=l...n. 
satisfies  the  nonl  inear  Langevin  equation 


<^(x,t)  =  fjO^x.tbx.t]  +  n^x.t) 


(2) 


The  dimension  of  space  is  d  and  the  dynamic  exponent  is  denoted  by  z.  Also, 


A(kw;kV)  -  A^kwikV)  -  - 

6^(ku))^(M 


(3) 


Cl  '  yg ' 

The  correlation  function  G(koj)  and  response  function  R(ku>)  are  given  by 
the  functional  derivatives: 


2 

^*-11 

6ijG(ku))  =  6h.(-k,-u)6h  .(k.oi)  'h,g  =  0  ’ 

'  J 

2 

(5  ^  W 

6ijR(kuj)  =  6h1(-k,-ajr5«ij(Vu)lh,g  =  0 


where  the  generating  functional  W  is  defined  by 


(4a) 


(4b) 


expW({h,g})  =  /[dij;1  ]/fd^21 


exp{-/dx/dt[iL(^1  (x,t) ,  i^(x,t),  ^(x,t)) 
-h( x,  t) (x,t)-g(x,t)^2(x,t)] 


(5a) 


with  L  =  l_i  +  L2  and 


L-j  =  ^2(x,t)  [ii»i  (x,t) 


f(^(x,t)  ,x,t)  -  ^  Z 


3f.  [ip-j(x,t),x,t] 


(5b) 


— 
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L2  =  -  J  Z  /dx  dx'/dt  dt '  Sjk(x,t;x',t,)^(x,t)v2(x,,t‘) 

J  >  k 

+  Ij-  z  /dx  dx'  dx-'/dt  dt'  dt"  S . k ( x , t ; x '  , t 1 ; x “ , t " ) 

i » . i- 

vj(x,t)^(x',t')^(x",t") 

+  ...  (5c) 


are  the  n-point  connected  Green  functions  of  n . . 


We  restrict  our  discussion  to  the  special  situation  of  a  time-dependent- 
Ginsburg-Landau  (TDGL)  model  with 


fi  =  -r(x»t)— f 


<SH 


(6) 


6v-|(x,t) 


< i  ( x » t ) ^ ^  ( x '  ,t‘  )>  =  2r(x,t)6(x-x‘  )6(t-t'  )6.  ■ 

Then,  the  action  integral  is: 
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where  the  n^^'s  are  combinatorial  factors,  and 


r^o  =  i(r+k  +icj/rk)  , 


r2l  =  i ( r+k  -iw/rk)  , 


r22  '  1/2rk  ‘ 


We  write  A  =  A^  +  A',  with  the  Gaussian  fixed  point 

AG  =  , d^^TT^l 2^1  ^ k -k » -a)) 

(2  it) 

+  r21  (k,ai)i|I2(k,w)-J1(  -k,-u>)  +  r22(k)02(k,w)-^2(-k,-w)]  ,  ( 

[From  here  on,  we  shall  explicitly  display  the  inner  products,  indices,  and 
summation  signs.] 

Associated  with  the  Gaussian  fixed  point  are  the  exponents 

4  if  r  =  vkd 

\  -  0  an<1  ZG  E  z  *  <2,  if  r  =  r.  ’  l 

k 

Then,  defining 

r^(k»“)  =  ^(k»w)[^(2-d)  -  z  6a2  -  k?k  -  z  w~],  ( 

the  linearization  of  (1)  about  Afi  reads, 

3 A '  _  ,  r  die  rdto  l  ?.  „i , ,  N  SA 1 

jr-  -  (d+z)A  +  / - -r/p-  T,  i  r  (koo) — - 

3e  (2TT)d  ^  a=l  i  =  l  a  6^(M 


+  tr  [Aj,-1  (q,w;-q,-u))A’ (q,u);-q, -(.))] 

(  2tt) 
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We  denote  the  RHS  of  (13)  by  L  and  define  eigenfunctions  0  with  eigenvalues 


LO  -  y<\ 


If  the  last  term  of  (13)  was  absent,  the  solutions  to  (14)  would  follow  immediately. 
Thus,  invoking  the  identity 

eABe~A  =  B  f  [A,B]  +  1[A,(A,B)]  +  —  (15) 

for  operators  A  and  B,  we  consider  a  transformation  of  (14)  generated  by  an 


operator  y 


(eyte  y)(eyC)  =  L'C  =  y;v  =  yey'\ 


such  that  the  second  functional  derivative  terms  are  eliminated.  Such  a  trans¬ 
formation  is  obtained  by  choosing  y  to  satisfy  the  commutator  relations, 

[y, tr  [Ar  1  (q,m;-Q,-(o)A(g,w;-q.-.01]  =  0,  ( 1  7a ' 

( 2  it  ) d  G 


r  -  (j  k  /-dtu  .,1  /  .  \  -» 

[y,;-— d.  or  >->•  I  ik,ui) — -J. - -] 

( 2uj )  '  mi  6iji^(k,m) 


tr  [Ar']  (q,u>;-q,~a).\(q,u>;-q, -«.>)], 
4  (2.:r)d  dr  b 


where  is  the  matrix  with  entires  .  ■  :_  ( q  -q ,  -  and 


-q ,  -'■■) 


,'''1,1t(q,:a)6i(j“(-q,-a) 
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For  y  of  the  form. 


—  -./j—  tr  [c(k ,w)A( k ,u»;-k, -w)  1 , 

(2.)d 


I  ]Q) 


where  c(k,.j)  is  the  matrix  with  entries  c^J(k,w)  =  6^.  .c  jP(  k ,\>) ,  the  relation  (17a) 
is  immediately  satisfied  for  any  choice  of  c^(,.  Thus,  c  is  completely  determined 
by  the  second  commutator  in  terms  of  the  partial  differential  equations: 


—  —  1  ^  ^  ^ 

(2+z+k9.+zu9  )c,,(k,u>)  =  -T- — — - ^  6 ( k - 1  ) , 

K  i>3  1  I  c.  i  /  _  \  i  L- 


f12(q,u>) 


(20a) 


—  i  t1-!  2 ( q  ,uj ) 

(2+k3k+zw9,o)c12(k,w)  =  jf  -  ~ 2 

i 


<5(k-l ) , 


( 20b  1 


1  r„ , ( q ,u) 

(2+k3.+Zay3  )c91(k,w)  =  y . . .  6(k-l), 

k  ^  J  2  !r12(q,W)r 


(20c) 


(2-z+k3k+Zu)3^)c22(k,w)  =  0. 


2nd) 


We  require  that  any  entry  c^  of  c  vanish  when  the  correspondi no  entry  of 

A  1  a  _  I 

( r-~  )  ,  =  (A _  )  ,(q,w;-q,-w)  vanishes.  Thus,  the  solution  to  (20)  may  be 

b  Up  b  Up 

represented  in  the  form. 


c(k,u)  =  y  0(l-k)(rT*1)(k,u)) 


(21  ) 


The  generalization  to  any  number  of  internal  fields  is  immediate, 
transformed  Eq.  (14)  thus  reads. 


The 


(drz)O'  +  MX  r1  (k,co) — I 


60' 


(2"T)d  2ir  ui  n  (S^(k,m) 


yo 


(2?) 


The  O'  so  determined  are  homogeneous  polynomials  in  v-  and  Label  tne  deur 

of  tne  f"  by  2n,  n  =  1,2,...  and  ^  by  m.  The  degree  of  0  in  v,  is  therefore 
2n-m.  The  O'  are  of  the  general  form 


d  k ,  ck»  i , 

1'  =  n  .  / - 

\  n  ( 2  ) 


d  k.,  d.w. 

,  2n  r  2n  , ,  ^ 

:  (2Tr)d'  2tt  v2n,m  k1  ’"1  ”  '  '  ’K2rV'2n' 


wnere  n 


4 1 1  k  1  ’“’l )  ‘  \z  (  k2  ’ ,d2 }  ‘  ‘ '  va2n _  1  (  k2n -  P"2n  - 1  ^ ^ (  k2n  ’  "2n } 


c  (  k,  +  .  .  .  +k.,r))  o(,1+...+m2n). 


is  a  normalization  constant  depending  on  (n,m)  and  the  pairing. 


I  It 

and  v.,  satisfies 

on 


[d+z+n(2-d) -iiiz-Z ( k  •  \  +Za-.h  )]v 


i  'k  •  Za  i  Vu  '  v2n  ,m  ^^ngn' 


Tig  v ~  's  are  tnus  homoneneous  polynomials  in  the  k.  and  the  of  decrees 

C  M  ,  ill  1  1 

r  and  s,  respecti ve ly ,  with  the  associated  eigenvalues 


y0’  =  2n*-d(l-n)  +  z(i-m)  -  r-s. 

J  2  n ,  m 


Slaving  developed  the  ei genoperators ,  we  can  now  expand  A  in  terms  of 


A  =  Ar  f  ;  V  a!/S  (i9r’s  . 
G  n,m  r,s  2n,m  2n,m 


ing  inner  product,  we  have 


2n,m  r  r ,  s  .  ,  ,  N  c  r.  s 

t----  -  !v9’  _  -  n-Sii - ( m-1  )<>z]a0’ 
d-  ?n,m  11  J  2n ,m 


y'  c  1  c  ’  \'  '  <’  ’ 

+  Vn.m  K’f'  a2n  *  ,m '  ,l2n  '  ,  rn '  1 ' 


VI !  -  ° 


where  on  and  6z  represent  corrections  to  the  respective  Gaussian  values  and  >n 
represents  the  nonlinear  terms  in  A1. 

*0  0 

We  no  suppose  in  d  =  4--  dimensions,  the  fixed  point  value  of  j  is  of 
order  :  while  all  other  a2^’^  are  of  higher  order  in  .  We  further  suppose  that 

,  9 

the  leading  corrections  to  on  and  cz  are  of  order  Thus, 


0,0  _  ,  0,0;0,0;0,0» 

4,1  '  "  a4 , 1  ;4 , 1  ;4 , 1 


(28a) 


*r,s  _  c2  r,s;0,0;0,0  ,  0,0;0,0;0,0, 

3 2 n , m  "  a2n,m;4,l ;4 ,1  a4, 1  ;4, 1 ;4,1  b 


( 28b) 


where  a^  =  0^(0 .  ,0^)  and  o  is  the  quadratic  part  of  P. 

2 

We  obtain  the  expressions  for  6n  and  6z  to  order  >  in  terms  of  a^^: 
( i )  Conserved  order  parameter 


(n=l ,  m=l ,  r=2,  s=0) 

xn  _  2,0;0,0;0,0/  0,0;0,0;0,0,~^2 

'  a2,l  ;4,1  ;41  U4,1;4,1;A,P  c 


r(n+2)/2(n+8)2]c2, 


(2qa) 


(n=l ,  m=2,  r=-2,  s=0) 


x _  _  -2 , 0 ; 0 , 0 ;0 , 0 /  0 , 0 ;0 ,0 ;0 ,0 , 

“  2,2;4,1  ;4 , 1  ^a4,l ;4,1  ;4,1  ’ 


-  5n  =  -6n 


(29b) 


( i i )  Non-conserved  order  parameter 


(n=l ,  m=l ,  r=2,  s=0) 

=  [(n+2)/2(n+8)2]c2,  (30a) 
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( n= 1 ,  m=2,  r=0 ,  s=0) 

6z  =  [6  tn(4/3)-l ]6n .  ( 30b 

C .  REMARKS 

We  have  demonstrated  that  the  problem  of  large-scale  stochast. i c i  ty  of  a 
classical  nonlinear  dynamical  system  can  be  discussed  conveniently  in  terms  of 
a  path  integral  formalism.  To  understand  the  full  stochastic  behavior  of  the 
system  u*. },  it  is  necessary  to  introduce  a  set  of  conjugate  momentum  field 
variables  (y9l.  The  relevant  correlation  and  response  functions  G  and  P  are 
obtained  as  functional  derivatives  of  the  generating  functional  W  [Eos.  4  and  •' 

For  classical  problems  in  critical  dynamics,  a  one-part icle-irreduc if .Y 
differential  renormal  ization-group  generator  is  introduced  [Eg.  (il\  >  is 
equation  can  be  solved  conveniently  in  terms  of  a  complete  set  of  Gaussian 
operators.  The  readers  are  urged  to  consult  the  Physical  Review  article  by 
Vvedensky,  Chang,  and  Nicoll  (1983)  for  additional  details. 

Large-scale  stochast  i  city  in  Hamitonian  systems  with  mu!  t  i -resonant  ♦»*»*•••. 
(e.g.  ,  *ne  problem  oc  stochastic  n eating  of  a  magnetized  plasma  by  la  rue  .r  ; 
coherent,  waves)  can  be  studied  in  terns  of  Poincare  naps.  The  threshold  c  an 

scale  s toe  ha ti city  has  been  shown  recently  to  be  amenable  to  a  renormal i za ' ion 

scrieme  (Escande.  198?)  akin  to  that  is  discussed  in  VII  A  and  B. 
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VIII.  FIELD  SWELLING  INSTABILITY  IN  ANISOTROPIC  PLASMAS  (BASU  AND  COPPI  1 Q8B ) 

A.  INTRODUCTION 

In  plasma  regimes  where  the  degree  of  col  1 isional i ty  is  relatively  low, 
the  particle  distribution  in  velocity  space  can  depart  considerably  from  a 
Maxwellian.  Then  plasma  collective  modes  can  be  excited  and,  as  a  result, 
induce  a  substantial  change  in  the  particle  distribution  in  velocity  space. 

Here  we  consider  the  case  where  the  electron  distribution  is  characterized  by 
an  anisotropic  temperature.  The  plasma  is  imbedded  in  a  relatively  weak  magnetic 
field  such  that  the  particle  pressure  is  of  the  same  order  as  the  maanetic  Dressure 
and  the  electron  temperature  anisotropy  (TgjL  f  T  )  is  relative  to  the  direction 
of  the  magnetic  field.  We  consider  collective  modes  that  can  be  described  by 
moment  (fluid-like)  equations,  and  so  do  not  rely  on  wave-particle  resonance 
processes  which  depend  on  the  evolution  of  a  small  portion  of  the  particle 
distribution  in  velocity  space.  The  plasma  motion  is  allowed  to  be  decoupled 
from  the  magnetic  field  lines.  That  is,  the  so-called  "frozen-in  condition" 
is  not  imposed.  Thus,  a  new  kind  of  instability  is  found  to  be  excited  when 
Tgi_  T  •  this  case  the  magnetic  field  is  perturbed  from  its  equilbrium 
state,  and  transverse  electron  thermal  energy  is  transported  toward  the  regions 
where  the  magnetic  field  is  weakened,  both  by  a  particle  flow  and  an  effective 
thermal  conductivity  along  the  magnetic  field  lines.  The  local  increase  in  the 
particle  transverse  pressure  tends  to  make  the  magnetic  field  "swell"  further 
locally  and  the  excited  mode  is  amplified  as  a  consequence.  We  point  out  that, 
as  shown  in  a  parallel  paper\  modes  of  the  type  treated  here  can  be  important 
in  producing  so-called  magnetic  reconnection  when  a  plasma  is  confined  in  a 
magnetic  field  configuration  that  contains  a  ("neutral")  surface  where  the  field 


vanishes.  Magnetic  reconnection  corresponds  to  a  change  of  the  field  topology, 

as  in  the  case  where  magnetic  "islands"  are  formed  out  of  a  configuration  that 

initially  has  straiaht  and  parallel  field  lines.  In  fact,  the  formation  of 

these  "islands"  would  not  occur  if  the  field  lines  were  constrained  to  move 

together  with  the  plasma.  The  analysis  of  neutral  sheet  configurations  where 

magnetic  reconnection  takes  place  is  of  special  interest  in  space  physics. 

We  note  that  fluid-like  instabilities  driven  by  a  plasma  temperature 

anisotropy  had  been  found  earlier,  and  the  one  which  is  closest  to  that  considered 

2 

in  this  paper  is  usually  referred  to  as  the  "mirror"  instability.  However,  this 
is  derived  under  the  assumptions  that:  a)  the  "frozen-in  condition"  is  valid  and 
b)  the  effective  thermal  conductivity  along  the  magnetic  field  is  negligible. 

In  particular,  the  threshold  value  of  Tei/Ten  above  which  this  instability  can 
be  excited  is  six  times  that  for  the  instability  we  describe  in  the  present 
paper. 


B .  A\AL_YSI_S 

We  write  the  equilibrium  electron  distribution  as 

T_.  1/2 


_  _  _n _ ,  eV'-  n(  Vi 

^ p  -  a/2  (t  )  exP(-  p-r 

(2"Ta  ,/.nJ3/Z  Te»  2Te. 


m  v 
e  •• 


FT 


)  , 


e^  e 


e" 


and  the  uniform  equilibrium  magnetic  field  as  =  Boz.  We  analyze  perturbations 

% 

of  the  form  B  •*  exp(-icot  +  i k Xx  +  ik„z)  with  ^  •  Q  and 


v thi  7k„I  ’  vthe  ’ 

:  eB  / m  c,  .J.  =  eB  /m  e,  v 


o  e 


i 


o  i 


the 


(2Te../me)1/?’  and  vthi  - 


We  introduce  the  vector  and  scalar  potentials  by  B  -  V  x  (A),  and  E  =  ( i.j/c )  A- . 

'V 

and  consider  the  quasi-neutrality  condition  n;  =  rv  f i . e . ,  ’'-A  -  pi. 


We  assume 


VIII-'? 


for  simplicity  that  the  ion  temperature  is  negligible  in  comparison  to  the  elec¬ 
tron  temperature.  Thus,  the  longitudinal  (cold)  ion  momentum  conservation 
equation  yields 


%  ie  £ 

u .  =  -  t 

i "  wm.j  z 


(1) 


From  the  transverse  (cold)  ion  momentum  conservation  equation  we  obtain,  for 


OJ  <  il .  , 


cM 


J  ~  m  A  +  JiL_  A  -  SL  __ 

uix  B  y  B  SI.  x  si.  B 
o  0  1  10 


(2) 


uiy 


ickJ  +  .i*L_  A 

Bo  BoSIi  y 


Bo  x 


Then  the  ion  mass  conservation  equation  gives 


ik.  ^  iek  ^ 

n  '  B  y  _  2  lz 

r\  J  — 


m^co 


(3) 


(4) 


'V/  f\  ^ 

for  u.  ituA  /B  that  corresponds  to  | k„/kA  |  >  |w/sl.  |.  In  order  to  proceed  to 
i  x  y  o  1 

% 

evaluate  ng  we  write  the  electron  pressure  tensor  as 


P 

-e 


P0  I  +  (P 


e„ 


PeJb  b  , 


where  I_  is  the  unit  dyadic,  b  =  B/B,  and  the  pe's  are  the  scalar  components. 

The  relevant  linearized  longitudinal  electron  momentum  conservation  equation  is 


n  = 


en 


'X, 

E. 


ik„ 


e,. 


+  ik..  (p„  -  PQ  ) 


% 

B. 


e„ 
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b  v  ien£z  A 

pe„  =  V.  "e  =  TT'  *  1k-  (pe„  '  ped  6„  > 


as  Tg  =0  given  that  w/k.,  ;  v^)e,  and 


n  p  ip  ^  i  k  p 

n  k"V,  z  Bo  Te„  * 

Then  the  quasi-neutrality  condition  gives 

wc2  eE„  Te  k  A 

(  1  -  — § — ) -I—  =  _^k  \  * 

u  2  '  k„T  T  B  y 
m  e„  e„  o 
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where  w  =  k„  T  /in..  From  the  transverse  electron  momentum  conservation 
s  e  M  i 

equation  we  find 


u  -  J-  A  (T  -  T  )  A 

ex  Bo  y  eB  2  e"  ex  x 
o 


Then,  from  the  Ampere's  law 


j  -  a 

x  4r  x 


we  obtain 


y  j  -  o  (B„  -  6J]  A 

k^v  d  * 

K  VA 


2  ck 
- - {  -  -i--i 

,2  2  '  *  >  ’ 

k  vA 


where  k2  =  k^  +  kA,  v?  =  /  (4unm.),  B„  =  Hrp  /B  2  and  ;-;A  =  rt-'P  / . 

M  0  I  t  u  v)  0 

Next,  we  consider 

ck'’  *'- 

J  =  r  A 
y  4tt  y 

and  derive  J  from  the  total  momentum  conservation  equation 


-  lion  m .  u  ■ 

1  i  x 


(k*  pex  ■  “•  (pc„  -  De,’  B*  ‘  c  Vo  • 
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Combining  Eqs.  (2),  (9)  and  (10)  we  obtain 

,2  k2  'v  4tt  k  % 

D  -  ~Y2 - T  -  6j]  A  =  i-y-^  p 

k  v„  2  k  y  k  B 


We  derive  p  from  the  linearized  equation  of  state 


ii*2peiv.^-peiv.,;en  =  ?.(i^;)  , 


where  V„  (BQ-V)/Bo  and 


1  m  /  dv  v?  (v„  f  -  u„  f  )  , 


4  r-v  III 

e„  2  e 


e  e„  e' 


f  being  the  perturbed  electron  distribution. 

% 

Since  we  have  not  found  an  easy  way  to  derive  q  x  from  a  'fluid1  approach, 

m 

we  calculate  it  directly  from  Eq.  (13)  after  solving  the  perturbed  linearized 
Vlasov  equation.  Introducing  the  polar  coordinates  0,v„  and  va  in  velocity 
space  we  obtain  the  moment 

'V  ^ 

2tf  ^  ?rrei  Ez  w-(l-oi  )k  V  E  v  k  V 

r  do  f  =  — —  t-JL  -  _ § — r _£  -  —k.  jL._i  a  q 

Jo  00  e  T  lk„  u)-k„v„  ‘k„  2c  fl  V 


v  k  v^  ^ 

•  l'-"e'«ylfe  •  "4> 

e 

ftp  ft, 

that  is  relevant  to  the  evaluation  of  u  and  q  x  .  Here  n  =  T  /T  and  we  find 

e„  e„  e  ej.  e„ 

'V  ft.. 

ft,  Tp  k  A  ck  B 

=  •  '"Te.  c;  if  -i"Te1(Te,.  '  V  e^V  '  (15) 
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O.  %  'Xj  V, 

Noting  that  2p  V-u  -p  V„uo  =  iwp  (n  /nj+ik^p  u  ,  we  obtain,  using  Eqs.  (5), 

“j.  c  cA  till  c  ca 

(7)  and  (15)  in  (12), 


Me  iP  -v-  e  N 

k"Te„  *  '  Te„  B„ 


and,  since  p  =  T  n  +  n  T  , 

ex  ex  e 


Te  B 

=  ( 1  T~^  JT  ’ 

ex  e„  o 


fhen  combining  Eqs.  (6),  (11)  and  (16)  we  finally  arrive  at  the  dispersion 


relation 


(:%  -  1)  [1 
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.  (2q  )  - 

2  2  2  T  ' >’x> 

k  vA  2k  e„ 
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'  2  v  r 
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r<  k2  Te 

»*>3*  -t  r;  - 

•v  c  (I 


It  should  be  mentioned  that  the  above  dispersion  relation  can  be  obtained  from 
tne  linearized  Vlasov-Maxwel 1  equations  by  evaluating  n  from  Eq.  (14)  and  u  w 


0  *'  fe  *  T-  ru-( 

6x  e 
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We  first  consider  k„/k  <<2/0„ .  Then  one  of  the  roots  of  Eq.  (18)  corre- 

2  2 

sponds  to  be  slow  magnetosonic  mode  with  w  'mj  ,  and  is  given  by 


\  [1  +  ^  (2  -  ^)]  -  1  +  Bx  (1  -  y*)  . 


This  becomes  unstable  if 


'e  8 

2  2  2 

The  other  root,  corresponding  to  the  fast  magnetosonic  mode  with  w  n,  k  v^,  is 


given  by 


2  R  Ip 

1  -  ~y~2  +  -f  (2  -  T^)  *  0  , 
k  e  ii 


and  become  unstable  if 
^e  B 


This  we  call  the  field-swelling  instability. 
In  the  other  limit,  k„»kJL,  we  have 


J  =  k^  Vfl  [1  -  \  (S„  -  0J]  , 


the  instability  criterion  being 


B„  -  Bx  >  2  , 


2  % 

and  this  is  the  so-called  'firehose'  instability.  In  this  case,  E_,  =  0, 

n.  %  n,  % 


B  =  0  implying  n  =  n.  =  0  and  T  =  0. 
z  3  e  i  ej. 
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Another  well  known  instability  driven  by  an  electron  temperature  anisotropy 
3 

is  the  Weibel  mode.  This  differs  from  the  one  we  present  here  in  that  it  develop 
only  in  zero  or  weak  magnetic  fields,  for  which  -  °  .  it  does  not  'pvnlve  *-nc 
ion  population,  it  does  not  produce  density  perturbations,  and  is  cnaracteri zed 
by  k  •  E  =  0.  Another  related  kind  of  instability  reported  in  the  literature 
involves  resonant  interaction  of  the  mode  having  u>  =  k„v^  and  k„>>kx  with  the 
tail  of  the  ion  distribution  in  velocity  space.  However,  in  this  case  the  driv¬ 
ing  temperature  anisotropy  is  that  of  the  ion  population. 

Finally,  we  refer  to  Eq.  (19)  and  note  that  for  T  =  T  and  >  1  we  have 

6  j.  £  ii 

2  2  2  2 

u)  t  -v'(2/B)  =  k„  v^.  Thus,  .a  ■*  0  as  Bq  ■+  0,  and  we  may  argue  that  this  is  the 

limit  from  which  the  familiar  "tearing"  mode,  that  is  driven  by  the  (spatial) 

gradient  of  the  current  density,  emerges  in  a  plane  "neutral"  sheet  confinuration. 

In  fact  the  influence  of  an  electron  temperature  anisotropy,  adding  the  effects 

of  the  Weibel  instability  to  the  theory  of  the  tearing  mode,  was  analyzed  in 

Ref.  5.  In  Ref.  1  it  is  pointed  out  that,  for  T  f  T  the  root  of  Eq.  (lol 

ex  e.. 

)  2 

■i  (1-T  /T  )/[l-T  / ( 2T  )]  does  not  vanisn  as  B  >0  ard  a  new  kind  of 

s  ex  e,i  cx  e  ii  o 

"fast"  magnetic  reconnection  process  associated  with  this  mode  can  take  place 
in  a  neutral  sheet  configuration. 

C .  CONCLUSION 

In  summary,  we  have  shown  that  collisionless  plasmas  imbedded  in  a  maonenc 
field  and  with  the  electron  temperature  transverse  to  the  magnetic  field  lamer 
than  the  longitudinal  temperature  can  be  subject  to  an  instability  tnat  produie- 
transport  of  transverse  electron  thermal  energy  into  the  regions  where  the  maaneti 
•iriu  weakened  by  the  perturbation  while  the  plasma  motion  is  decoupled  from 
that  of  the  magnetic  field  lines. 
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